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ABSTRACT.  We  show  how  estimates  on  the  Hausdorff  dimensions  of  the  Holder 
singularities  of  a  function  can  be  derived  from  the  properties  of  its  expansion 
on  an  orthonormal  wavelet  basis.  The  validity  of  the  multifractal  formalism 
is  examined  and  results  concerning  the  genericity  of  multifractality  in  a  given 
function  space  are  deduced. 


Multifractal  analysis  is  a  recent  field,  introduced  in  the  context  of  turbulence  in 
the  mid  80s;  its  purpose  is  to  analyze  the  pointwise  Holder  regularity  of  functions, 
and  to  understand  how  this  regularity  fluctuates  from  point  to  point.  One  tries 
to  determine  the  spectrum  of  singularities  of  the  function  /  studied:  It  is  a  novel 
function  which  associates  to  each  Holder  exponent  H  the  dimension  of  the  sets  of 
points  where  /  has  this  given  pointwise  regularity. 

We  will  give  a  detailed  account  of  the  interaction  between  wavelet  analysis  and 
multifractal  analysis.  This  interaction  was  made  possible  because  of  two  remarkable 
properties  of  wavelet  expansions: 

•  It  is  possible  to  characterize  the  pointwise  Holder  regularity  of  a  function 
by  simple  estimates  on  the  decay  rate  of  its  wavelet  coefficients,  corre¬ 
sponding  to  the  wavelets  localized  near  the  point  considered. 

•  Wavelets  are  unconditional  bases  of  “most”  function  spaces. 

The  second  point  is  relevant  because  a  central  problem  in  multifractal  analysis 
is  to  relate  the  spectrum  of  singularities  of  a  function  /  to  the  function  spaces  which 
contain  /.  Formulas  which  perform  this  bridge  are  called  multifractal  formalisms. 

This  paper  is  split  into  five  parts: 

Section  1  is  introductory;  we  define  the  pointwise  Holder  regularity  and  the 
different  notions  of  dimension  that  we  will  use.  These  notions  are  illustrated  by  the 
study  of  two  functions:  The  Weierstrass  ‘“nowhere”  differentiable  functions,  and 
the  devil’s  staircase.  These  toy  examples  allow  us  to  introduce,  in  a  very  simple 
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setting,  some  of  the  tools  that  will  be  explored  and  used  in  the  following:  The 
wavelet  criterion  for  irregularity  and  the  mass  distribution  principle. 

Section  2  is  devoted  to  the  wavelet  characterization  of  the  pointwise  Holder 
exponent  and  the  relationship  between  Holder  regularity  and  local  oscillation.  We 
give  a  new  formulation  of  this  criterion  in  terms  of  local  suprema  of  the  wavelet 
coefficients,  the  wavelet  leaders,  which  paves  the  way  for  the  new  multifractal 
formalism  introduced  in  the  following  section. 

Section  3  is  devoted  to  the  study  of  this  new  multifractal  formalism.  First, 
a  heuristic  justification  is  given;  though  it  follows  the  lines  of  the  initial  argument 
given  by  Parisi  and  Frisch,  this  formulation  is  given  in  terms  of  wavelet  leaders  and 
not,  as  usual,  in  terms  of  either  increments,  or  wavelet  coefficients.  This  has  the 
advantage  of  eliminating  some  causes  of  failure  of  the  multifractal  formalism.  We 
prove  that  this  multifractal  formalism  yields,  for  all  functions  that  have  a  minimal 
uniform  regularity,  an  upper  bound  for  their  spectrum  of  singularities;  we  give  an 
alternative  derivation  of  the  multifractal  formalism  in  terms  of  local  oscillations, 
and  we  prove  that  both  methods  are  equivalent. 

Section  4  is  a  complement  to  Section  3:  We  derive  properties  satisfied  by  the 
scaling  function  which  appear  in  the  formulation  of  the  multifractal  formalism.  We 
show  that  it  indicates  which  oscillation  spaces  contain  the  function  considered  and 
we  give  some  properties  of  these  spaces. 

Section  5  gives  “generic”  results  of  multifractality.  First,  we  discuss  what 
genericity  can  mean  in  infinite  dimensional  function  spaces.  We  focus  on  one  partic¬ 
ularly  important  notion:  Prevalence.  Finally,  we  perform  the  multifractal  analysis 
of  “almost  every”  function  (in  the  sense  of  prevalence)  belonging  to  a  given  Besov 
or  oscillation  space. 


1.  What  is  multifractal  analysis? 

In  this  first  section  we  introduce  the  main  concepts  that  are  relevant  in  mul¬ 
tifractal  analysis,  and  we  perform  the  “multifractal  analysis”  of  two  very  simple 
examples:  the  Weierstrass  functions  and  the  devil’s  staircase. 

1.1.  Pointwise  regularity.  Multifractal  functions  are  used  as  models  for  sig¬ 
nals  whose  regularity  may  change  abruptly  from  one  point  to  the  next.  Our  first 
task  is  therefore  to  define  what  is  meant  by  pointwise  regularity.  It  is  a  way  to 
quantify,  using  a  positive  parameter  a,  the  fact  that  the  graph  of  a  function  may 
be  more  or  less  “ruguous”  at  a  point  xo- 

Definition  1.  Let  a  be  a  nonnegative  real  number,  and  xo  E  Md;  a  function 
/  :  — >■  M  is  Ca(xo)  if  there  exists  C  >  0,  S  >  0  and  a  polynomial  P  of  degree  at 
most  [a]  such  that 

(1.1)  if  \x  —  xo\  <6,  \f(x)—P(x  —  xo)\<C\x  —  xo\a- 

Remarks:  We  use  the  (unusual)  convention  [a]  =  a  —  1  if  a  is  an  integer.  This 
convention  implies  that  the  polynomial  P  is  unique  even  when  a  is  an  integer;  the 
constant  term  of  P  is  necessarily  f(x o);  P  is  called  the  Taylor  expansion  of  /  at  xo 
of  order  a. 
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We  will  also  use  a  slightly  weaker  notion;  /  is  C^og{x o)  if  there  exists  C  >  0, 
d  >  0  and  a  polynomial  P  of  degree  at  most  [a]  such  that 

(1.2)  if  \x  —  xo\  <  S,  | f(x)  —  P(x  —  *o)|  <  C \x  —  x0\a  log  (l/\x  —  *o|)  • 
Definition  2.  The  Holder  exponent  of  /  at  xo  is 

hf(xo)  =  sup  {a  :  /  is  C'a(x o)}  . 

Note  that  (1.1)  implies  that  /  is  bounded  in  a  neighbourhood  of  xo;  therefore, 
the  Holder  exponent  is  defined  only  for  locally  bounded  functions.  The  Holder 
exponent  is  defined  point  by  point  and  describes  the  local  regularity  variations  of 
/.  Some  functions  have  a  constant  Holder  exponent;  they  are  called  monoholder 
functions.  We  will  see  that  it  is  the  case  of  the  Weierstrass  functions;  it  is  also 
the  case  of  the  sample  paths  of  the  Brownian  motion  which  satisfy  almost  surely: 
Vx,  hsix)  =  1/2.  Such  functions  display  a  “very  regular  irregularity”.  As  a  first 
example,  we  now  study  the  pointwise  regularity  of  the  Weierstrass  functions;  indeed, 
this  will  allow  us  to  see,  in  a  very  simple  setting,  how  wavelet  techniques  can  be 
used  to  prove  irregularity. 

1.2.  The  Weierstrass  functions.  The  Weierstrass  functions  are  defined  by 

OO 

Wa,b(x)  =  cos(5n*)> 

n  =  1 

where  B  is  assumed  to  be  larger  than  1,  so  that  the  series  is  lacunary,  and  A  is 
assumed  to  be  smaller  than  1,  so  that  the  series  converges  normally.  Derivating 
term  by  term,  it  is  clear  that,  if  AB  <  1,  Wa,b  est  differentiable.  We  assume  that 
AB  >  1  which  will  actually  imply  that  Wa,b  is  nowhere  differentiable. 

Proposition  1.  If  A  <  1  <  AB,  the  Holder  exponent  of  Wa,b  is  a  constant 
function  which  is  equal  to  a  =  —  log  A/log  B. 

Proof  of  Proposition  1:  Let  us  first  check  that  Wa,b  is  C'a(x o)  for  any 
xo  E  M.  In  the  difference 

OO 

Wa,b(x)  ~  Wa,b{x o)  =  Y,AU  (cos (Bnx)  -  cos {Bnx0))  , 

n  —  1 

we  can  either  bound  the  difference  of  cosines  simply  by  2  or,  using  the  mean  value 
theorem,  by  Bn\x  —  xo  |.  Let 

N  =  r-iog(k-Eoi)' 

log  B 

Using  the  first  bound  for  n  >  N  and  the  second  one  for  n  <  N ,  we  get 

N  oo 

\Wa,b(x)  -  Wa,b(xo)\ <  +  2  An . 

n—l  n=N 

We  have  to  sum  up  two  geometric  series.  Because  of  the  value  taken  for  N ,  the 
first  sum  is  bounded  by  C(AB)N  <  C\x  —  xo\a ,  and  the  second  one  is  bounded  by 
CAN  <  C\x  —  xo\a ,  hence  the  regularity  result  holds. 

In  order  to  prove  irregularity  of  Wa,b  at  every  point,  we  will  use  a  technique 
which  predates  the  wavelet  methods  developed  in  Section  2. 
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Definition  3.  Suppose  that  ip(x)  satisfies  the  following  conditions 

C 


C  f 

IV’(*)I  <  "j - j — 12  et  /  ip{x)dx  =  0. 

1  +  \x\  J  K 


Let  /  G  L°°(M),  a  >  0  and  b  G  M;  the  continuous  wavelet  transform  of  /  is  dehned 

by 


(1.3) 


C{a,b)  =  f  f(x)ip 

J  K 


x  —  b\  dx 


Lemma  1.  If  /  G  C'a(x o)  for  a  G  [0, 1],  then  the  continuous  wavelet  transform 
of  /  satisfies 

\C{a,b)\<C{aa  +  \b-x0\a). 

Proof  of  Lemma  1:  Since  ip  has  a  vanishing  integral, 

' x  —  b\  dx 


that 


C(a,b)  =  f  (f(x)  -  f(x0))ip 

J  K 


\C(a,  b)  |  <  C  [  \x-x0\a  ip  ( 

J  K  V 


x  —  b 


dx 

a 


<  C 


\x  —  b\a  +  \b  —  x0\a  dx 


1  + 


v-b  I 


<C{aa  +  \b-x0\a). 

Let  us  come  back  to  Proposition  1.  The  idea  of  its  proof  is  to  pick  ip  so  that  the 
computation  of  C'(a,  b)  may  be  as  simple  as  possible.  It  is  the  case  if,  in  (1.3),  only 
one  frequency  of  the  Weierstrass  function  is  selected.  A  possible  choice  is  therefore 
to  pick  for  ip  a  function  whose  Fourier  transform  ip  is  C 2  and  satisfies 


supp(ip)  C 


et  ^(1)  =  2. 


The  hypotheses  of  Lemma  1  are  clearly  satisfied  and,  if  a  =  B  N 


C(B~N ,  b) 


Vd"  cos  (Bnx)iP 

n  =  1 


x  —  b 


B~N 


dx 


B~N 


=  Vd"  /  cos (Bn~Nu  +  Bnb)iP(u)du. 

n  =  1 

Because  of  our  choice  for  ip,  the  only  nonvanishig  integral  corresponds  to  n  =  N  so 
that 

C(B-N,b)  =  ANetBNb; 

thus, 

if  a  =  B~n  ,  then  \C(a,b)\  =  aa. 

Lemma  1  implies  that  the  Holder  exponent  cannot  be  larger  than  a  at  any  point 
and  Proposition  1  follows.  In  particular,  the  Weierstrass  functions  are  nowhere 
differentiable. 
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1.3.  Measure  and  Hausdorff  dimension.  Multifractal  analysis  studies  func¬ 
tions  whose  Holder  exponent  can  jump  from  one  point  to  the  next.  If  such  is  the 
case,  the  numerical  computation  of  the  Holder  exponent  of  a  signal  is  completely 
instable,  and  also  quite  meaningless.  One  rather  wishes  to  obtain  a  more  global 
information  such  as:  Does  the  Holder  exponent  take  a  given  value  HI  And,  if  such 
is  the  case,  what  is  the  size  of  the  set  of  points  where  hj  takes  this  value?  This 
question  rises  the  following  problem:  Which  kind  of  ‘size’  should  we  use?  The  an¬ 
swer  can  be  justified  by  the  study  of  many  mathematical  functions:  Usually,  there 
is  a  ‘most  probable’  exponent,  which  is  taken  almost  everywhere;  therefore,  ‘size’ 
cannot  mean  ‘Lebesgue  measure’,  which  would  not  allow  to  distinguish  between 
the  sizes  of  all  but  one  sets  of  points  where  a  given  Holder  exponent  appears.  The 
Hausdorff  dimension  is  more  fitted  for  this  purpose.  Let  us  first  recall  the  notion 
of  d-dimensional  Hausdorff  measure,  see  [33,  34]. 

Definition  4.  Let  A  c  Md.  If  e  >  0  and  d  £  [0,  d\,  we  denote 


where  R  is  an  e-covering  of  A,  i.e.  a  covering  of  A  by  bounded  sets  of 

diameters  \Ai  \  <  e.  The  inhmum  is  therefore  taken  on  all  e-coverings. 

For  any  d  £  [0,  d\,  the  d-dimensional  Hausdorff  measure  of  A  is 

mess  (A)  =  lim  Ml5. 

'  '  e-MI  e 

(Note  that  the  limit  exists  (it  can  take  the  value  +oo)  since  M /  is  a  decreasing 
function  of  e.) 

Definition  5.  Let  A  C  Md;  then,  there  exists  do  £  [0,  d]  such  that 

Vd  <  do,  mess  (A)  =  +oo 
Vd  >  do,  mess  (A)  =  0. 

This  critical  do  is  called  the  Hausdorff  dimension  of  A. 

The  existence  of  do  is  a  consequence  of  two  straightforward  assertions:  If  d  <  S' , 

mess  (A)  <  +oo  =y  mess'{A)  =  0  and  mess'{A)  >  0  =>■  mess(A)  =  +oo. 

In  order  to  bound  the  Hausdorff  dimension  of  a  set,  it  is  sufficient  to  consider 
particular  coverings.  By  contrast,  it  is  harder  to  obtain  a  lower  bound  directly 
from  the  definition,  since  it  requires  to  consider  all  possible  coverings.  The  fol¬ 
lowing  “mass  distribution  principle”  replaces  the  study  of  all  e-coverings  by  the 
construction  of  one  measure  (which  is  often  naturally  supplied  by  the  fractal  set 
considered). 

Proposition  2.  Let  m  be  a  probability  measure  supported  by  A  C  Md.  As¬ 
sume  that  there  exist  d  £  [0,rf],  C  >  0  and  e  >  0  such  that,  for  any  ball  B  of 
diameter  less  than  e, 

m{B)  <  C\B\S. 

Then  mess  (A)  >  1/C. 

Proof:  Let  {Bi}  be  an  arbitrary  e-covering  of  A.  Then 

1  =  m{A)  =  m  (IJB,-)  <  ]Tm(B0  <  \B{\S . 
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Other  definitions  of  fractal  dimension  have  been  introduced  and  are  used  in 
multifractal  analysis.  The  simplest  notion  is  supplied  by  the  box  dimensions. 


Definition  6.  Let  A  C  Md;  if  e  >  0,  let  Ne(A)  be  the  smallest  number  of  sets 
of  radius  e  required  to  cover  A. 

The  upper  box  dimension  of  A  is 


dims  (A)  =  limsup 

e-MI 


log  Ne(A) 

-loge 


The  lower  box  dimension  of  A  is 


dimnSA)  =  liminf 

-  '  '  e-MI 


log  Ne(A) 

-loge 


One  important  drawback  of  the  box  dimensions  is  that,  if  A  is  dense,  then  the 
box  dimensions  take  invariably  the  value  d.  Since  most  multifractal  functions  of 
interest  have  dense  sets  of  Holder  singularities,  box  dimensions  are  unable  to  draw 
any  distinction  between  the  sizes  of  these  sets.  This  explains  why  box  dimensions 
are  not  used  in  the  definition  of  the  spectrum  of  singularities. 

Another  notion  of  dimension  which  retains  some  flavour  of  the  box  dimension, 
but  has  better  mathematical  properties,  is  supplied  by  the  packing  dimension. 
Invented  by  C.  Tricot,  see  [99],  it  can  be  defined  by 


dim 


■P  (-4)  =  inf 


1  sup  (  dim} 


■A,  :  A  C  \J  Ai 


8  =  1 


(the  inhmum  is  taken  over  all  possible  partitions  of  A  into  a  countable  collection 
Ai).  The  packing  dimension  has  been  extensively  used  in  the  context  of  multi¬ 
fractal  measures;  however  its  use  is  much  more  scarce  for  functions;  see,  however, 
Proposition  8  where  an  upper  bound  for  packing  dimensions  of  Holder  singularities 
is  obtained;  see  also  [50]  where  other  bounds  of  packing  dimensions  are  obtained 
for  a  different  notion  of  Holder  singularity. 


1.4.  The  spectrum  of  singularities.  The  Holder  exponent  of  multifractal 
functions  may  take  a  given  value  on  a  fractal  set.  If  such  is  the  case,  one  wishes  to 
determine  the  Hausdorff  dimension  of  this  set. 

Definition  7.  Let  /  :  and  H  >  0.  If  if  is  a  value  taken  by  the 

function  hf(x),  let 

(1.4)  EH  =  {x0£Rd:  hf{x0)  =  H}. 

The  spectrum  of  singularities  (or  Holder  spectrum)  of  /  is 

df(H)  =  dim(Eji) 

(we  use  the  convention  df(H)  =  — oo  if  H  is  not  a  Holder  exponent  of  /).  The 
support  of  the  spectrum  is  the  set  of  values  of  H  for  which  Eh  0. 

Remarks: 

•  We  will  also  consider  the  set 

Gh  =  {x0£  Rd:  f^Cfog(x „)}. 


(1.5) 
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•  These  global  notions  can  also  be  defined  locally;  if  SI  C  is  a  nonempty 
open  set,  we  define 

Eg  =  Eg  fl  SI  and  d^(H)  =  dim(Eg); 
clearly,  dj(H)  =  sup  d^(H). 

•  The  spectrum  of  singularities  is  defined  on  M  +  U{+  oo}. 

If  hf  takes  at  least  two  finite  values,  /  is  said  to  be  mutifractal.  If  hj  takes  only 
one  finite  value,  /  is  said  to  be  monofractal ;  one  example  is  the  devil’s  staircase, 
considered  in  Section  1.5.  In  the  examples  we  will  consider,  df(H)  will  usually 
take  nonnegative  values  on  a  whole  interval  [Hmin,  Hmax\.  If  Hmin  yS  Hmax,  its 
computation  requires  the  study  of  an  infinite  number  of  fractal  sets  Eg]  in  such 
cases,  the  term  ‘multifractal’  is  completely  justified. 

To  perform  the  multifractal  analysis  of  a  function  means  to  determine 
its  spectrum  of  singularities. 

One  can  meet  two  types  of  multifractal  functions.  A  first  one  is  supplied  by 
inhomogeneous  functions,  which  are  smoother  in  some  regions  than  in  others.  This 
case  is  often  met  in  image  analysis.  Indeed,  an  image  is  a  patchwork  of  textures 
with  different  characteristics.  Their  spectrum  of  singularities  reflects  the  mono-  (or 
multi-)  fractal  nature  of  each  component,  and  also  of  the  boundaries  (which  may 
also  be  fractal)  where  discontinuities  appear.  In  such  situations,  the  determination 
of  the  local  spectrum  of  singularities  d^(H)  for  the  different  ‘components’  SI  is 
more  relevant.  By  contrast,  homogeneous  multifractal  functions  present  the  same 
characteristics  everywhere.  The  following  definition  makes  this  notion  precise. 


Definition  8.  Let  /  be  a  multifractal  function  (therefore,  its  Holder  exponent 
takes  at  least  two  values);  /  is  an  homogeneous  multifractal  function  if  the  spectrum 
of  singularities  of  its  restriction  to  any  nonempty  open  set  SI  is  independent  of  SI 
(and  therefore  coincides  with  the  whole  spectrum  of  singularities  of  /). 

1.5.  The  devil’s  staircase.  The  devil’s  staircase  is  probably  the  simplest 
fucntion  whose  multifractal  analysis  can  be  completed  in  an  elementary  way.  Its 
Holder  singularities  are  located  on  a  fractal  set,  the  triadic  Cantor  set,  which  is 
defined  as  follows. 

Let  x  E  [0,  1];  x  can  be  written  (in  base  3)  as 

OO 

x  =  J2~^j  with  aj-e  {0,1,2}. 
j= i 

The  triadic  Cantor  set  K  is  the  set  of  x  such  that  ay  £  {0,  2}  for  all  j. 

The  devil’s  staircase  is  defined  as  follows:  If  x  K,  at  least  one  of  the  ay  is 
equal  to  1.  Let  l  =  inf{j  :  ay  =  1};  then 


i-i 


v(x)  =  y 


3  =  1 


27 +  1 


a d_ 
2' ' 


The  function  T>(x)  is  thus  defined  almost  everywhere  on  [0,  1].  It  is  then  extended 
by  continuity  on  [0,  1]. 

One  easily  checks  that  V  is  increasing  and  that,  on  K,  its  Holder  exponent 
takes  the  value  log  2/ log  3;  it  is  equal  to  +oo  at  the  other  points. 
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In  order  to  compute  the  spectrum  of  singularities  of  V ,  we  have  to  determine 
the  Hausdorff  dimension  of  K.  Let  e  >  0  and  n  be  such  that  3_n  <  e  <  3_n+1. 
Since  K  can  be  embedded  in  the  union  of  2n  intervals  of  length  3_n,  using  these 
intervals  as  e-covering,  we  obtain  that  the  Hausdorff  dimension  of  K  is  at  most  d. 

A  lower  bound  for  the  Hausdorff  dimension  of  K  is  obtained  through  an  appli¬ 
cation  of  the  mass  distribution  principle.  Since  V  is  increasing  from  [0,  1]  to  [0,  1], 
its  derivative  fi  is  a  probability  measure.  Furthermore,  since  V  is  locally  constant 
on  the  complement  of  K,  then  Supp(pi)  =  K.  Let  I  =  [xo,x]  C  [0,  1].  The  Hoder 
regularity  of  V  implies  that 

p(I)  =  V(x)  -  V(x0)  <  2|/|los2/log3. 

Therefore,  msp(K)  >  and  the  mass  distribution  principle  yields  that 


dim(K)  > 


log  2 
log  3' 


Hence,  we  have  proved  the  following  statement. 


Proposition  3.  The  devil’s  staircase  is  monofractal;  its  spectrum  of  singular¬ 
ities  is  supported  on  ^  {+oo}  where 


dj) 


log  2 
log  3 


,  ^  and  d(+oo)  =  1. 
log  3 


1.6.  Notes.  The  notion  of  pointwise  Holder  regularity  has  been  generalized 
in  several  ways.  A  first  one  was  introduced  by  A.  P.  Calderon  and  A.  Zygmund  in 
1961,  and  leads  to  a  weaker  condition,  see  [23]:  /  belongs  to  Tp(x o)  if  there  exists 
a  polynomial  P(x  —  x o)  of  degree  at  most  [a]  such  that 

(1.6)  ( ~7  [  \f{x)  -  P(x  -  x0)\p  dx\  <Cra. 

\r  JB(x0,r)  J 

Clearly,  if  /  £  Ca(x o),  then  Vp  >  1,  /  E  T?(x o).  This  definition  is  particularly 
useful  when  one  deals  with  functions  which  are  not  locally  bounded.  Indeed,  (1.6) 
makes  sense  as  soon  as  /  E  Lploc',  therefore,  it  is  a  natural  substitute  for  pointwise 
Holder  regularity  when  functions  in  Lploc  are  considered.  This  is  particularly  rele¬ 
vant  when  dealing  with  applications  where  the  signal  studied  is  not  locally  bounded; 
for  instance,  it  is  the  case  in  fully  developed  turbulence  where  singularities  of  Holder 
exponent  -1  corresponding  to  thin  vorticity  filaments  can  be  observed,  and  in  mam¬ 
mography  images,  where  microcalcihcations  also  appear  as  singularities  of  Holder 
exponent  -1,  see  [5].  If  /  E  Lpoc,  one  defines  the  p-exponent  at  xo  by 

hf,P(x 0)  =  sup{a  :  feTp(x0)j, 
and  the  p-spectrum  of  /  by 

df,P(H)  =  dim({xo  '■  hftP(x0)  =  H}). 

These  notions  were  studied  by  C.  Melot  in  [85]. 

Let  us  mention  another  situation  where  the  p-exponent  is  more  relevant  than 
the  Holder  exponent.  Suppose  that  S2  is  a  domain  with  a  fractal  boundary;  up  to 
now,  only  one  parameter  was  available  to  classify  these  domains:  The  dimension 
of  their  boundary.  At  each  point  xo  of  the  boundary,  the  Holder  exponent  of  the 
characteristic  function  /  =  In  takes  the  value  0;  this  is  in  sharp  contrast  with  the 
1-exponent  which  can  clearly  take  any  value  in  [0,  oo],  depending  on  the  behavior  of 
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the  boundary  near  xo',  this  remark  allows  one  to  perform  a  multifractal  analysis  of 
fractal  boundaries  and  leads  to  a  whole  function  (the  1-spectrum)  as  discriminating 
parameter  between  these  domains;  this  yields  a  a  much  richer  classification  tool  for 
fractal  domains,  see  [59]- 

On  the  other  hand,  Holder  regularity  can  be  strengthened  as  follows,  see  [50]. 
We  denote  by  meas(A)  the  Lebesgue  measure  of  a  set  A.  Let  a  >  — d;  a  point  xo 
is  a  strong  a-smgularity  of  /  if  there  exist  C,  C"  >  0  such  that  V  P  polynomial  of 
degree  at  most  a,  Vj,  3 Aj,  Bj 

(  meas(Aj)  >  C2_dj,  meas(Bj)  >  C2_dj 

<  Vr  6  Aj  U  Bj,  \x  —  xo\  <  2_J 

{  Vt  e  Aj,  Vj/  e  Bj,  (f(x)  -  P(x  -  xo ))  -  (f(y)  -  P(y  -  x0))  >  C' 2~aP 

Note  that,  if  a  <  1,  the  last  condition  reduces  to  f(x)  —  f(y)  >  C"2~aj .  Bounds 
on  the  packing  dimension  of  the  strong  a-singularities  have  been  obtained  in  [50]; 
see  also  [59]  where  estimates  on  the  the  packing  dimension  of  the  boundary  of  a 
domain  in  M.d  are  obtained  in  terms  of  the  wavelet  coefficients  of  the  characteristic 
function  of  the  domain. 

A  natural  question  is  to  determine  which  nonnegative  functions  h(x)  are  Holder 
exponents.  An  exact  characterization  of  the  Holder  exponents  of  continuous  func¬ 
tions  was  obtained  by  P.  Anderson,  see  [2]  (following  results  of  S.  Jaffard  [49]  and, 
independently  of  K.  Daoudi,  J.  Levy-Vehel  and  Y.  Meyer  [30]):  h(x)  is  the  Holder 
exponent  of  a  continuous  function  f  if  and  only  if  it  is  the  lim  inf  of  a  seguence  of 
continuous  functions.  (Note  that  the  problem  is  still  open  if  /  is  not  assumed  to  be 
continuous) .  The  first  constructive  proofs  used  deterministic  constructions  for  /. 
In  practice,  for  simulation  purposes,  one  needs  to  construct  random  processes.  This 
led  to  the  introduction  of  mult ifr actional  Brownian  motions,  which  are  Gaussian 
processes  generalizing  the  Brownian  motion.  This  study  was  initiated  by  A.  Be- 
nassi,  S.  Jaffard  et  D.  Roux  in  [18],  and  independently  by  J.  Levy-Vehel  et  R.  F. 
Peltier  in  [76].  The  most  general  constructions  were  obtained  by  A.  Ayache  and  J. 
Levy-Vehel,  see  [10].  One  can  also  consult  the  survey  paper  [16]  by  A.  Benassi,  S. 
Cohen  and  J.  Istas,  which  covers  a  wider  range  of  related  problems. 

The  notion  of  Holder  exponent  is  adapted  to  functions  whose  regularity  changes 
abruptly  from  point  to  point.  When  it  is  not  the  case,  the  more  stable  notion  of 
local  Holder  exponent  can  be  used: 

Hf(x 0)  =  inf{a  :  3d  >  0,  /  £  Ca([x0  -  d,  x0  +  d])}. 

J.  Levy-Vehel  has  studied  the  properties  of  this  exponent.  In  particular,  in  collab¬ 
oration  with  S.  Seuret,  he  proved  the  following  characterization:  The  local  Holder 
exponents  of  continuous  functions  are  exactly  the  nonnegative,  lower  semicontin- 
uous  functions,  see  [77].  However,  up  to  now,  this  notion  has  had  no  impact  on 
multifractal  analysis:  Indeed,  all  multifractal  functions  we  will  meet  have  a  constant 
local  Holder  exponent. 

The  idea  of  associating  fractals  to  measures  or  functions  can  be  traced  back 
to  the  books  by  Mandelbrot  in  the  70s  and  80s  [82,  83]  ;  see  also  [84]  where  the 
main  contributions  of  B.  Mandelbrot  to  multifractal  analysis  are  collected,  and  [32] 
where  C.  Evertsz  and  B.  Mandelbrot  give  an  extremely  pedagogical  introduction 
to  multifractal  measures.  Multifractals  were  introduced  by  G.  Parisi  and  U.  Frisch 
in  [92]  for  modeling  fully  developed  turbulence,  following  the  pioneering  work  of  B. 
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Mandelbrot;  indeed,  B.  Mandelbrot  had  introduced  multiplicative  cascades  mod¬ 
els  in  [81],  which  were  studied  by  J.-P.  Kahane  et  J.  Peyriere  in  [66]  (see  also 
[14,  15,  89]  and  references  therein);  it  is  remarkable  that  these  cascades,  initially 
introduced  as  models  for  the  dissipation  of  energy  in  turbulent  flows  before  the  no¬ 
tion  of  multifractal  was  introduced  actually  supplied  the  first  mathematical  models 
of  multifractal  measures. 

Here  is  a  very  simple  example  of  a  multiplicative  cascade  which  is  a  multifractal 
measure.  It  is  constructed  recursively  on  [0,1]  as  follows.  Let  a  £  (0,1),  b  = 
l  —  a,  and  p  £  (0, 1).  We  pick  at  random  /u([0,  1/2])  =  a  or  b,  with  probabilities 
respectively  p  and  1  —  p,  and  p((  1/2,  1])  =  1  —  /u([0,  1/2]);  Once  the  measure  of 
a  dyadic  interval  A  has  been  determined,  we  split  this  interval  into  two  “sons” 
of  the  same  length;  the  measure  of  its  left  son  A'  is  also  picked  at  random  and 
takes  the  value  p( A')  =  ap( A)  or  (1  —  a)p(X)  with  probabilities  respectively  p 
and  1  —  p,  and  the  measure  of  its  right  son  is  p(X)  —  p(X').  One  thus  constructs 
iteratively  a  random  probability  measure  on  [0,  1].  This  is  the  simplest  example 
of  a  multiplicative  random  cascade;  there  exists  a  huge  literature  on  multifractal 
cascades,  see  [14,  35,  89]  and  references  therein  for  a  mathematical  analysis  of 
these  cascades,  and  [101,  12,  15]  for  important  recent  generalizations. 

We  won’t  consider  multifractal  measures  in  these  notes;  let  us  just  mention  that 
the  notion  of  Holder  exponent  is  replaced  by  the  local  dimension  of  the  measure  p 
defined  as 

/>„(*„)  =  i,mmflog("l[,r/,;,°+li])). 

<5->-0  log  (d) 

see  [19,  21,  32,  44,  91,  34],  However,  a  relationship  between  Holder  regularity 
for  measures  and  functions  can  be  explicitely  established  in  the  one-dimensional 
setting:  If  p  is  a  probability  measure  supported  on  M,  let  f(x)  =  p((—o o,  *]);  one 
immediately  checks  that 

if  hll{x0)  £  [0,  1),  then  hfi(x0)  =  hf(x 0). 

The  problem  we  mentioned  concerning  the  Holder  exponent  can  also  be  raised 
for  the  spectrum  of  singularities:  Which  functions  d(H)  are  spectra  of  singularities 
of  continuous  functions?  We  only  have  partial  answers:  limsup  of  nonnegative 
piecewise  constant  functions  are  spectra,  see  [50],  but  we  do  not  know  if  this  class 
covers  all  possible  spectra.  Furthermore,  one  would  rather  wish  to  characterize  the 
spectra  of  homogeneous  multifractal  functions;  in  this  case,  only  very  partial  results 
are  known,  see  [55]. 

There  exists  many  generalizations  of  Cantor’s  triadic  set  and  of  the  devil’s 
staircase;  see  for  instance  [65]  where  more  general  perfect  sets  are  constructed  and 
their  relationships  with  trigonometric  series  are  studied.  Another  generalization, 
fractal  strings,  was  introduced  by  M.  Lapidus  and  C.  Pomerance  in  [71,  72];  it 
was  then  used  in  a  number  of  papers,  including  those  by  M.  Lapidus  and  H.  Maier 
[69,  70],  C.  He  and  M.  Lapidus,  see  [42,  43]  and  M.  Lapidus  and  C.  Pomerance, 
see  [73]  ;  see  also  the  reference  book  [74]  by  M.  Lapidus  and  M.  van  Frankenhuysen. 

The  box  dimensions  yield  an  important  information  on  the  fractal  geometry  of 
the  graph  of  functions,  which  is  independent  of  the  spectrum  of  singularities.  Gen¬ 
eral  formulas  which  allow  to  derive  the  box  dimensions  of  the  graphs  of  arbitrary 
functions  from  their  wavelet  expansions  are  available,  see  [53]  and  [67].  The  prob¬ 
lem  of  determining  the  Hausdorff  dimensions  of  graphs  is  often  very  hard,  especially 
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in  deterministic  settings;  see  however  [17],  and  [11,  95,  96]  where  estimates  are 
derived  from  wavelet  expansions. 

The  notion  of  homogeneous  multifractal  functions  is  important  for  modelling; 
indeed,  with  regards  to  fully  developed  turbulence,  Parisi  and  Frisch  conjectured 
that  its  behavior  is  universal,  i.e.  that  the  spectrum  of  singularities  does  not  de¬ 
pend  on  the  particular  fluid  considered,  on  the  limit  conditions,...  In  particular,  if 
this  hypothesis  is  verified,  the  spectrum  of  singularities  of  the  velocity  of  a  turbu¬ 
lent  fluid  should  be  independent  of  the  region  where  the  fluid  is  being  inspected, 
and  therefore,  the  velocity  of  a  turbulent  fluid  is  expected  to  be  an  homogeneous 
multifractal  function.  Usually,  the  mathematical  functions  which  are  known  to  be 
multifractal  are  actually  homogeneous  multifractals.  We  will  see  several  examples 
backing  this  assertion. 

2.  Wavelets  and  Holder  regularity 

Orthonormal  wavelet  bases  are  a  privileged  tool  to  study  multifractal  functions 
for  several  reasons.  A  first  one,  exposed  in  this  section,  is  that  the  Holder  exponent 
can  be  characterized  by  simple  local  decay  conditions  on  the  wavelet  coefficients. 
Another  reason,  discussed  in  Section  3,  concerns  the  study  of  the  multifractal 
formalism.  It  is  a  formula  which  is  used  to  derive  the  spectrum  of  singularities  of 
a  function  /  from  the  knowledge  of  some  function  spaces  which  contain  /.  Here 
again,  wavelets  play  a  key  role  since  these  spaces  are  defined  by  simple  conditions 
on  the  wavelet  coefficients. 

We  will  just  recall  some  properties  of  orthonormal  wavelet  bases  that  will  be 
useful  in  the  following.  We  refer  the  reader  to  [27,  31,  64,  75,  80,  87]  for  detailed 
expositions  of  this  subject. 

2.1.  Orthonormal  wavelet  bases.  Orthonormal  wavelet  bases  are  construc¬ 
ted  through  the  help  of  a  multiresolution  analysis.  It  is  a  sequence  of  closed  sub¬ 
spaces  of  T2(Md)  denoted  by  Vj  ( j  £  Z)  and  satisfying: 

(1)  Vj  £  Z,  Vj  c  vj+1, 

(2)  Vj  £  Z,  f(x)  £  Vj  «=>  /( 2x)  £  Vj+1, 

(3)  3<p(x)  £  Vo  such  that  the  functions  tp(x-k)  (k  £  Zd)  form  an  orthonormal 
basis  of  Vo . 

(4)  n  u  =  {0}  and  (^J  Vj  is  dense  in  L2(Rd). 

The  multiresolution  analysis  is  r-smooth  if  ip  is  C'r  and  if  the  daip,  for  |a|  <  r,  have 
fast  decay. 

The  subspace  Wj  is  defined  as  the  orthogonal  complement  of  Vj  in  Vj+i-  Under 
these  assumptions,  there  exist  2d  —  1  functions  satisfying  the  same  regularity 
and  decay  properties  as  ip  and  such  that  the  (x  —  k)  (i  =  1,  ...,  2d  —  1,  k  £  Zd) 
form  an  orthonormal  basis  of  Wo-  Using  Point  2,  it  follows  that  Vj  £  Z,  the 
2dV’1f)V)  (23 x  —  k)  (k  £  Vd)  form  an  orthonormal  basis  of  Wj. 

By  construction,  the  Wj  are  orthogonal,  and  we  can  use  two  possible  decom¬ 
positions  of  T2(Md)  as  a  direct  orthogonal  sum 

+  oo  +oo 

L2  =  0  Wj  or  L2  =  V 3  +  0  Wj, 

3=0 


3  =  ~  °° 
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which  leads  to  two  possible  ways  to  write  a  function  of  L 2 

(2.1)  /w  =  EEE4/')(2Ji-*). 

j£%k£hd  i 

or 

oo 

(2.2)  f(x)  =  Ck<p{x-k)  +  YY  Y0),^  ^  x  ~ 

k£hd  j=0 k£hd  i 

the  Cj  k  are  the  wavelet  coefficients  of  / 

(2.3)  c)  k  =  2^  f  f{x)^(2jx  -  k)dx, 

J  Kd 

and 

(2.4)  Ck  =  f  f(x)<p(x-k)dx. 

JWd 

A  consequence  of  Point  2  is  that  the  functions  2d^2p(2^  x  —  k)  form  an  orthonormal 
basis  of  Vj .  Thus,  if  Pj(f)  denotes  the  orthogonal  projection  of  /  on  Vj, 

Pj  (f)(x)  =  Yj  Cjikip(22x  —  k)  where  Cjtk  =  2_dj  f  f{x)ip(22x  —  k)dx. 

k  ,d 

Remarks:  In  (2.1)  or  (2.2),  we  do  not  choose  the  L 2  normalisation  for  the 
wavelets,  but  rather  an  L°°  normalisation  which  is  better  fitted  to  the  study  of 
Holder  regularity. 

Note  that  (2.3)  and  (2.4)  make  sense  even  if  /  does  not  belong  to  L 2;  indeed, 
if  one  uses  smooth  enough  wavelets,  these  formulas  can  be  interpreted  as  a  duality 
product  betweeen  smooth  functions  (the  wavelets)  and  distributions  (/).  We  will 
see  the  cases  of  Sobolev  and  Besov  spaces  in  Section  4.1. 

A  simple  example  of  multiresolution  analysis  is  obtained,  in  dimension  1,  by 
taking  $  =  l[o,i]!  in  this  case,  Vj  is  the  space  of  square  integrable  functions  which 
are  constant  on  each  interval  [^-,  and  a  possible  choice  for  ’k  is  1  [o,  1/2]  —  l[i/2,i]  • 

The  corresponding  wavelet  basis  is  called  the  Ptaar  basts. 

Among  the  families  of  wavelet  bases  that  exist,  two  will  be  particularly  useful 
for  us: 

•  Lemarie-Meyer  wavelets,  such  that  ip  and  both  belong  to  the  Schwartz 
class; 

•  Daubechies  wavelets,  such  that  the  functions  ip  and  can  be  chosen 
arbitrarily  smooth  and  with  compact  support. 

We  will  often  assume  that  the  multiresolution  analysis  is  of  tensor  product  type 
i.e. ,  if  x  =  (xi,  .  .  .Xd),  then 

ip(x)  =  $(x1)$(x2)  ■  ■  .®(xd), 

where  $(*)  is  associated  to  a  1-dimensional  multiresolution  analysis.  If  ^(a:) 
denotes  the  corresponding  1-dimensional  wavelet,  we  can  take  for  d-dimensional 
wavelets  the  functions 

//*)(*)  =  ^,l(*l)^,2(*2)  •  •  .^d(*d), 

where  denotes  either  $  or  'P,  and  where  the  choice  $(*i)$(*2)  •  •  •<f’(*d)  is  the 
only  one  excluded  (thus  there  are  indeed  2d  —  1  wavelets). 
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If  the  multiresolution  analysis  is  r-smooth,  the  wavelets  have  a  corresponding 
number  of  vanishing  moments,  see  [87]: 

If  | cl |  <  r,  then  f  ipl'l\x)xadx  =  0. 

JRd 

Therefore,  if  the  wavelets  are  in  the  Schwartz  class,  all  their  moments  vanish. 

2.2.  Wavelets  and  uniform  regularity.  We  start  by  characterizing  uniform 
regularity  properties  in  terms  of  the  wavelet  coefficients.  Let  us  recall  the  definition 
of  the  Holder  Ca  spaces. 

Definition  9.  Let  a  £  (0, 1).  A  function  /  belongs  to  C"(Md)  if 
/  £  L°°  (Md)  and  if  3C  >  0  such  that 

(2.5)  Vx,y£Rd  \f(x)-f(y)\  <C\x-y\a. 

Let  a  >  1  and  a  ^  N.  A  function  /  belongs  to  C"(Md)  if  f  £  L°°  (Md)  and  has 
partial  derivatives  of  any  order  \f3\  <  [a]  and  if,  when  \f3\  =  [a],  then 

(2.6)  Vx,y£Rd  \dPf(x)  -d^f(y)  \  <C\x  -y\a~^. 

The  Ca  norm  of  /  is  the  sum  of  ||  /  ((oo  and  of  the  inhmum  of  the  Cs  such 
that  (2.5)  or  (2.6)  holds.  The  following  proposition  shows  that  the  space  Ca  is 
characterized  by  a  very  simple  decay  condition  on  the  wavelet  coefficients. 

Proposition  4.  Let  a  >  0  be  such  that  a  LJ  and  assume  that  the  multires¬ 
olution  analysis  is  of  tensor  product  type  and  r-smooth  with  r  >  [a]  +  1.  Then 

f  \/k£Zd  \ck\<c 

(2.7)  /  e  Ca{Rd)  <^>300:  l 

{  Vj>0,Vi,  \/k£Zd  \c)>k\ <  C2~aK 

Remarks:  In  the  following,  we  will  usually  not  mention  the  regularity  needed 
for  the  wavelets,  which  will  be  assumed  to  be  “smooth  enough”;  the  minimal  reg¬ 
ularity  required  following  easily  from  the  computations. 

Proposition  4  holds  for  any  wavelet  basis  which  is  r-smooth  with  r  >  a  because 
the  matrix  of  a  change  of  wavelet  basis  is  continuous  on  the  space  of  coefficients 
satisfying  (2.7),  see  [87]. 


Proof:  Suppose  that  /  £  C"(Md);  then, 

\Ck\<  [  \f{x)\  \<p{x  -  k)\dx  <  C  ||  /  ||oo  • 

The  wavelet  ip ^  (x)  can  be  written  'Li(ti)'I,2(t2)  •  •  •  'Sd(xd),  where  the  if;  are  equal 
to  $  or  if,  and  at  least  one  of  them  is  equal  to  if.  Suppose  that  it  is  the  case  for 
’Ll.  If  k  =  (k1;  .  .  .,kd), 

cljk=2dj  /  /(*!,  .  .  . ,  Td)^(2JT1  -  ki)  .  .  ,^d(2Jxd  -  kd)dxi  .  .  ,dxd. 

The  function  iL(t)  has  its  first  [a]  +  l  moments  vanishing;  therefore  it  can  be  written 

*<*>  =  ThW 

where  if  still  has  fast  decay  and  a  vanishing  integral.  Let 

tp(x  1,  ■  ■  .Xd)  =  ^(ti)^2(t2)  •  •  -^d(xd)- 
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Integrating  by  parts  [a]  times  in  the  direction  xi,  we  get 

c)>k  =  2d3  J  f{x)rp{i\2ix  -  k)dx 

r  r/M  f 

=  (-1)M2*2-[“]j  /  — ±{x)ip(2 Jx-k)dx 
d  dx  2 

=  «-1»["l2<’2"["1'  /  (0M  -  (f))  ~  k)il 

(because  ip  has  a  vanishing  integral);  therefore,  using  the  decay  of  ip, 


c)k  \<  C 2dj2~^j 


2J 


a—  [a] 


c 


(1+  I  23 X  -  k  |)d+ 2 


dx  <  C2~a3. 


Let  us  now  prove  the  converse  result.  Assume  that 
(2.8)  \Ck\  <  C  and  |c*-jfc|  <  C2~aj . 

Let 

f-i(x)  =  £CW(*  -  k) 


and,  if  j  >  0, 

(2.9) 

let  us  consider  the  series 

(2.10) 


i  k 


+  oo 


£  fj 

j=- 1 


Using  (2.8)  and  the  decay  of  ip  and  ip O,  one  obtains  that  (2.9)  converges  uni¬ 
formly  on  any  compact  to  a  limit  fj  which  has  the  same  regularity  as  the  wavelets. 
Furthermore, 


2  ~a3 


Vj>  1,  \fj(x)  \<CJ2  (i+\2jx-k\)d+i 


<  C' 2~aj: 


therefore,  the  series  (2.10)  converges  normally  to  a  function  which  we  denote  mo¬ 
mentarily  by  g  and  which  belongs  to  L°° .  Similarly,  using  the  decay  of  the  d13  (p 
and  of  the  d13  ip1-1'1 ,  one  obtains 

(2.11)  if  /3\  <  [a]  +  1,  |  dpfj(x)  |<  U2^l-“)i. 

We  can  differentiate  term  by  term  the  series  fj  up  to  the  order  [a] ,  and  still 
obtain  a  normally  convergent  series;  thus  g  belongs  to  and 

V|/?|  <  [«],  I d/3g{x)\<C. 

Let  x  and  y  in  Rd  and  jo  be  defined  by  2_j0_1  <|  x  —  xq  \<  2 -j0;  if  \(]\  =  [a],  then 
I  X3 g(x)  —  d13 g(y)  |  is  bounded  by 

£  I  dPfjix)  -  d?f3(y)  |  +  £  |  dPfjix)  |  +  £  |  dPfjiy)  |  . 
j<j  0  j>j  0  j>j  0 
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Using  (2.11)  and  the  mean  value  theorem,  the  first  term  is  bounded  by 


£  I x~y 

3<3o 


sup  sup  \&>m 
|/3|=[a]  +  !  [x,x0] 


<C\x-y  | 

3  <3 o 


<C\x-x o  I 


Coming  back  to  the  definition  of  jo,  it  follows  that  the  first  term  is  bounded  by 
C  |  x  —  xq  |“.  The  second  and  third  terms  are  bounded  by 


<C\x-x o 

i>io 


a  — [a] 


J 

Let  us  now  prove  that  /  =  g.  We  consider  the  sequence  Fj  —  f  —  £  fj;  its 

j=- 1 

wavelet  coefhcients  vanish  for  j  <  J  and  it  shares  the  same  wavelet  coefficients  as 
/  if  j  >  J.  Therefore,  the  same  computation  as  above  shows  that  ||  Fj  ||oo<  C2~J , 
so  that  /  =  g;  hence  the  converse  result  holds. 


2.3.  Characterization  of  the  Holder  exponent.  We  start  by  introducing 
some  definitions  and  notations.  A  dyadic  cube  of  scale  j  is  a  cube  of  the  form 


k\  fci  +  U  kd  L  +  Q 

J  x  x  |_v,— iff- J  ’ 


where  k  =  (k\,  .  .  .kd)  E  ZC  Instead  of  indexing  the  wavelets  and  wavelet  coeffi¬ 
cients  with  the  three  indices  (i,j,  k),  we  will  use  dyadic  cubes.  Since  i  takes  2d  —  1 
values,  we  can  assume  that  it  takes  values  in  {0,  \}d  —  (0,  .  .  . ,  0);  we  will  use  the 
following  notations: 


• A  (=A(uu))  + 

•  c^  =  cik 

•  ip\(x)  =  UU  2-C  —  k), 


0, 


1 

T+I 


d 


k 

•  =  — . 

The  wavelet  ip\  is  essentially  localized  near  the  cube  A;  more  precisely,  when  the 
wavelets  are  compactly  supported 


3C  >  0  such  that  Vi,  j,  k,  supp  (ip\)  C  C  ■  A 

(where  C  ■  A  denotes  the  cube  of  same  center  as  A  and  C  times  wider).  Finally,  A  j 
will  denote  the  set  of  dyadic  cubes  A  which  index  a  wavelet  of  scale  j,  i.e.  wavelets 
of  the  form  tp\(x)  =  UH  2V  —  k)  (note  that  A  j  is  a  subset  of  the  dyadic  cubes  of 
side  2J+1).  We  take  for  norm  on  M.d 

if  x=(xi,...,  xd),  \x\=  sup  |  a?  *  | ; 


so  that  the  diameter  of  a  dyadic  cube  of  side  2  J  is  2  J . 


Definition  10.  The  wavelet  leaders  are  defined  by 
(2.12)  dx  =  sup  |cA'|- 

A'C  A 
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Note  that  if  /  E  L 00 ,  then 

|ca|  <  2*  J  \f{x)\\il>x{x)\dx  <  C sup \f{x)\; 
so  that  the  wavelet  leaders  are  finite. 

Definition  11.  Two  dyadic  cubes  Ai  et  A2  are  called  adjacent  if  they  are  at 
the  same  scale  and  if  dist( Ai,  A2)  =  0  (note  that  a  dyadic  cube  is  adjacent  to  itself). 
We  denote  by  Aj(xo)  the  dyadic  cube  of  side  2~J  containing  xo  and  by  Adj( A)  the 
set  of  3d  dyadic  cubes  adjacent  to  A.  Then 

dj(x  0)  =  sup  dy. 

\'£Adj(\j(x0)) 

Definition  12.  The  influence  cone  at  xo  consists  of  the  set  of  dyadic  cubes 
of  the  form  Aj(xo)  and  their  adjacent  cubes. 

The  wavelet  characterization  of  the  Holder  exponent  requires  the  following 
regularity  hypothesis,  which  is  slightly  stronger  than  continuity. 

Definition  13.  A  function  /  is  uniformly  Holder  if  there  exists  e  >  0  such 
that  /  E  Ce(Md). 

The  following  theorem  allows  to  characterize  the  pointwise  regularity  by  a  decay 
condition  of  the  dj(x 0)  when  j  — >■  +00. 

Theorem  1.  Let  a  >  0.  If  /  is  Ca(x 0),  then  there  exists  C  >  0  such  that 

(2.13)  Mj  >  0,  djixo)  <C2~aj. 

Conversely,  if  (2.13)  holds  and  if  /  is  uniformly  Holder,  then  /  belongs  to  C^og{x 0). 

Remark:  If  /  is  uniformly  Holder,  its  regularity  at  xo  is  therefore  determined 
by  the  behavior  of  the  d\  in  the  influence  cone  at  xo-  More  precisely,  one  can  notice 
that  the  characterization  supplied  by  Proposition  4  can  be  rewritten 

3 C  VA  dx  <  C2~aj. 

Therefore  pointwise  and  uniform  Holder  regularity  are  characterized  by  the  same 
decay  condition,  either  written  locally  (in  the  influence  cone)  or  uniformly. 


Proof  of  Theorem  1:  If  /  is  C'a(x 0),  then  there  exists  a  polynomial  P  of 
degree  at  most  [a]  such  that 

| f(x)  -  P(x  -  *o)|  <  C\x  -  x0\a. 

Let  j  >  0;  if  A'  C  3Aj(xo),  then,  using  the  vanishing  moments  of  the  wavelets, 
c\i  =  2dj  J  f(x)ip^'1  (2J  x  —  k’)dx  =  2dj  J (f(x)  —  P(x  —  xo))ip^'1  (2J  x  —  k')dx; 
thus  |ca'|  is  bounded  by 

C2dj'  f  \x-x0\a  <  C2dj<  [  I  X  -  n\'  |“  +  I  n\'  —  x0 


(1+  |  2J'x  -  k '  |)d+“+ 1 


(1+  |  2i'x  -  k'  \)d+a+ 1 


-dx 


(because,  Va,  a,  b  >  0,  (a  +  b)a  <  2"(a"  +  &"))•  The  change  of  variable  t  =  2J  x  —  k' 
yields 


(2.14) 


\cx,\<C(2~ai  +|/iA'  ~x0\a). 
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Since  j'  >  j  —  1  and  | fly  —  xo\  <  4d2  we  have  |ca'|  <  C 2  “•?,  so  that  dj(x o)  < 
C2~aP 


Let  us  now  prove  the  converse  result.  First  note  that,  since  /  is  uniformly 
Flolder,  Proposition  4  implies  that  the  wavelet  series  of  /  converges  uniformly  to  / 
on  any  compact.  Let  j'  be  given.  We  will  first  estimate  the  size  of  fji  and  of  its 
partial  derivatives.  If  A'  is  a  cube  of  side  2-J  ,  denote  by  A  (=  A(A'))  the  dyadic 
cube  defined  by 

•  If  A'  C  3Aj'(to),  then  A  =  Aji(xo), 

•  else,  if  j  =  sup{/  :  A'  C  3A;(*o)},  then  A  =  Aj(xo)  (and  it  follows  that 
2~J~1  <  \nx'  ~  *o|  <  4d2“J). 

In  the  first  case,  by  hypothesis,  |ca'|  <  dji(x o)  <  C. 2_"J  ,  and  the  sum  on  the 
corresponding  A'  satisfies  (as  in  the  proof  of  Proposition  4) 

<  C2~aj‘ . 

A' 

In  the  second  case,  |cv|  <  dj(x o)  <  C'2~aj  <  C\xo  —  H\'\a  and  the  sum  on  the 
corresponding  values  of  A'  of  scale  j'  satisfies 

<  fST'  \xo  ~ 

-  2-j  (1  +  1 2 j ' a?  -  k'\)d+a+1 

\x  —  xo\  -\-  \x  —  fi\i\ 

—  2—j  (x  _|_  |2 j'x  —  k'\)d+a+1 

<C(\x-x0\a  +2-ai'). 

Thus 

(2.15)  \fr(x)\<c{2-ad'  +  \x-x0\a); 

Similarly,  one  obtains,  for  any  [3  such  that  \f3\  is  smaller  than  the  regularity  of  the 
wavelets, 

(2.16)  | dpfj>(x)\  <  C2^j‘  (2~ajl  +  |a?  —  a?0|“)  • 

OO 

In  particular,  if  \j3\  <  [a],  the  series  fj(x  o)  converges  absolutely;  (recall  that, 

j=- 1 

by  convention,  if  a  is  an  integer,  [a]  =  a  —  1)  and  we  can  define  the  polynomial 

P(x-x0)=  J2  ^  £  dfifj{x0y, 

|/3|<[a]  ’  i  =  -l 

P  is  the  sum  of  the  polynomials 

Pj{x-x o)  =  ^  dl3fj(x o); 

|/3|<[a]  1 

Since  Pj  is  the  Taylor  expansion  of  fj  of  degree  [a]  at  xo , 

\fj(x)  -  Pj(x  -  x0)\ <  C\x  -  x0\[a]+1  sup  ^  \dP  fj\ 

[x,X0]  |/3|_[q,]  +  1 


^ca'V’a'(t) 

A' 
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(2.17)  <  C\x  -  *0|W+12([“]+1)i  ( 2~aj  +  |*  -  z0|“)  • 

We  will  now  bound 

(2.18)  \f(x)  -  P(x  -  x0)\=  ^2  (fj(x)  -  Pj(x  -  x 0)) 

j>- 1 

By  hypothesis,  there  exists  e  >  0  such  that  /  £  Ce(Md).  Let  jo  and  j i  be  defined 

by 

2“io-i  x  _  Xq  |<  2_-Jo  and  j\  =  — 1 — 

[  e 

We  have  jo  <  C|  log  |a?  —  a?o  1 1 -  If  j  <  jo,  using  (2.17),  the  sum  of  the  corresponding 
terms  of  (2.18)  is  bounded  by 

Y  C\x  -  xolM+^M+ib  (2-“f  +  I*  -  xo\a)  <  Cjo\x  -  *o|“. 

3  <3 o 

If  jo  <  j  <  ji,  (2.15)  implies  that  \fj(x)\  <  C*k  —  «o|“;  so  that 

J  1  j 

l/j (a?) |  <  C'ji  |a?  -  *o|“  <  —  k  -  *o|"  |log(k  -  *o|)|  • 

j=jo 

If  j  >  ji,  since  /  is  uniformly  Holder,  Proposition  4  implies  that  |jj(*)|  <  C 2_ej, 
and 

OO 

X]  |/j(z)|  <  C2-£f1  <  C2-“f°  <  Ck-^oT- 

i=ii 

Finally,  using  (2.16), 

^i^k-*o)i<  cy  Y 

3>3o  j>3o\/3\<[a] 

<CJ2  Y 

3>3 o  |/3|<[a] 

<  C\x  —  X o|". 

The  converse  part  of  Theorem  1  is  thus  proved. 

2.4.  Quasisure  Holder  regularity  in  C"(IBk)-  We  will  now  prove  a  quasi- 
sure  result  which  is  a  direct  consequence  of  Theorem  1.  Let  us  first  recall  Baire’s 
category  theorem: 

If  if  is  a  complete  metric  space,  every  countable  intersection  of  open 
dense  sets  is  dense. 

If  a  property  V(x)  holds  (at  least)  on  a  countable  intersection  of  open  dense 
sets,  it  is  said  to  hold  quasi-surely.  Note  that  a  quasi-sure  property  V  does  not 
necessarily  hold  on  a  large  set;  indeed,  in  ld  a  quasi-sure  property  may  hold  only 
on  a  set  of  measure  0  (and  even  of  dimension  0).  It  rather  means  that  the  property 
holds  on  a  “very  dense”  set.  Indeed,  countable  intersections  of  such  sets  are  still 
dense. 
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Proposition  5.  Let  a  >  0  and  such  that  a  N.  The  functions  of  C"(Md)  are 
quasi-surely  mono-Holder  with  Holder  exponent  h(x)  =  a.  (The  set  of  mono-Holder 
functions  of  exponent  a  even  contains  a  dense  open  set  of  C"(Md)). 

Proof  of  Proposition  5:  We  assume  that  a  smooth  enough  wavelet  basis  has 
been  chosen.  Any  function  /  £  Ca  can  be  written  in  the  form 

f(x)  =  ^2ck<p(x  -  k)  +  ~  k ), 

k  i  j>0  k 

where 

sup(|cfe|,  2"J  \clj  k\)  <  C. 
j,k 

We  can  choose  for  equivalent  norm  in  C"(Md)  the  inhmum  of  the  Cs  such  that  this 
inequality  holds.  If  TV  is  a  given  integer,  let  En  be  the  subset  of  C  "(Md)  dehned 
by:  F  £  En  if  all  \ck\  and  2"-J|c*  fe|  are  nonvanishing  integer  multiples  of  2~N . 
Clearly 

V/  £  C“(Md)  dist(f,  En)  <  2~N , 

so  that  U£=1  En  is  dense  in  In  a  metric  space  E,  we  denote  by  B(x,  R) 

the  open  ball  of  center  x  and  radius  R.  Let  Fjy  =  En  +  B( 0 ,  2— JV— 1 ) .  Fn  is  an 
open  set.  Any  function  /  of  Fn  satisfies 

Vj,*  2“J'|c*->fc|  >  2~n~1 

and  the  direct  part  of  Theorem  1  implies  that  the  Holder  exponent  of  /  is  everywhere 
at  most  a.  Since,  on  the  other  hand,  any  function  /  of  Fn  belongs  to  C“(M  d),  f  is 
clearly  a  mono-Holder  function  of  exponent  a.  The  union  of  the  Fn  is  the  expected 
open  dense  set. 

2.5.  Holder  exponent  and  local  oscillations.  We  will  now  see  a  conse¬ 
quence  of  Theorem  1  which  supplies  an  alternative  definition  for  Holder  regularity. 
In  this  new  definition,  the  polynomial  P ,  which  may  be  difficult  to  exhibit,  does 
not  appear  explicitly. 

Let  /  :  — y  M  and  h  £  Md.  The  first  order  difference  of  /  is 

{Ahf){x)  =  fix  +  h)  -  f{x). 

If  n  >  1,  the  differences  of  order  n  are  dehned  recursively  by 

(Anhf)(x)  =  (Anh-1f)(x  +  h)-(Anh~1f)(x). 

Definition  14.  If  A  is  a  convex  subset  of  Md,  the  oscillation  of  order  n  of/ 
on  A  is 

OSnf{A)=  sup  \{Alf){x)\. 

[x,x-\-nh\cA 

Note  that 

OS*  (A)  =  sup  f(x)  —  inf  f(x). 

x£A  X^A 

Proposition  6.  If  /  £  Ca(x 0),  then 
(2.19)  Vn  >  [a]  +  1  Ve  >  0  \OSnf{B{x0,e))\<Cea. 

Conversely,  if  /  is  uniformly  Holder  and  if  (2.19)  holds,  then  /  £  C“og(x o). 
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Note  that,  if  0  <  a  <  1,  the  triangular  inequality  implies  that 

f&Ca{x 0)  «=>  Ve>0  |OS}(5(*0,e))l  <  Cea. 

Proof  of  Proposition  6:  Assume  that  /  £  C'a(x o).  Since 
|  f(x)  -  P(x  -  *o)|  <  c\x  -  x0\a, 

it  follows  that 

if  \x  —  xo\  <  e,  then  |A h(f(x)  ~  P{%  —  *o))|  <  2nCe". 

But,  if  n  >  deg(P),  A^P(x  —  xo)  =  0.  It  follows  that  |A^/(*)|  <  Cea. 

In  order  to  prove  the  converse  result,  we  need  the  following  property  of  one¬ 
dimensional  wavelets. 

Lemma  2.  Assuming  that  the  multiresolution  analysis  is  n-smooth,  for  any 
n  >  1,  there  exists  a  function  6n  with  fast  decay  and  such  that  ip  =  A^28n. 
Furthermore,  if  ip  is  compactly  supported,  we  can  pick  6n  also  compactly  supported. 

Proof  of  Lemma  2:  Recall  that,  in  the  construction  of  one-dimensional 
wavelets  by  multiresolution  analysis,  the  embedding  V-i  C  Vo  implies  the  exis¬ 
tence  of  coefficients  4  such  that 

9  (|)  =J2lkv,(x  - 

k 

which  can  be  rewritten  in  the  Fourier  domain  as 

£(2£)  =  m0{p)Lp{p), 

where  mo  is  C°°  and  27r-periodic.  Furthermore,  mo  vanishes  at  n  to  the  order  at 
least  r,  see  [87].  One  can  pick  for  wavelet  ip  the  function 

f(20  =  +  7 r)£(£). 

It  follows  that  ip( 2£)  vanishes  at  the  multiples  of  27T  to  the  order  at  least  r.  Fur¬ 
thermore,  if  <p  is  compactly  supported,  mo  is  a  trigonometric  polynomial,  so  that 
(after  perhaps  translating  ip),  fho(£.  +  tt)  can  be  written 

(2.20)  m0(£  +  7r)  =  (e*-l)rP(O, 

where  P  is  a  trigonometric  polynomial,  and  therefore 

n 

p(0  =  J2  a^ki- 

k  =  —  n 


In  dimension  1,  to  apply  A)]  to  /  amounts  to  multiply  /(£)  by  elh ^  —  1.  Therefore, 
we  look  for  a  function  6n  such  that 

m  =  (e^l2-l  )"M0- 

We  can  dehne  6n  by 


MO 


-i£/2  mo(£  +  Q 

(gie/2  _ 


9 


Since  the  function  n(£) 


m0{£,  +  7r) 
(e*V/2  -  l)n 


is  C°° 


and  47r-periodic,  6n  has  fast  decay. 


Furthermore,  if  ip  is  compactly  supported,  (2.20)  implies  that  n(£)  is  a  trigonomet¬ 


ric  polynomial,  and  therefore  6n  is  compactly  supported.  Hence  Lemma  2  holds. 
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Let  us  now  prove  the  converse  part  in  Proposition  6.  Assume  that,  in  dimension 
d,  we  use  the  tensor-product  wavelet  basis.  Then 

UT)  =  ^(*1)  •  •  ■'&d{xd), 

where  the  MP  are  either  ip  or  <p,  but  at  least  one  of  them  is  equal  to  ip.  Suppose  for 
instance,  that  =  ip{x\).  Then 

c)tk  =  2dj  f  2-U  -  h)  . .  .U(2Vd  -  kd)dx 


=  j  ^,1(*i)  •  •  ■'$d{xd)dx. 

Clearly,  in  dimension  1, 

(2.21)  J  f(x)(Akg)(x)dx  =  j {^-hf){x)9{x)dx; 

so  that 


(2.22)  c)tk  =  {-lT  J  (Mei/2/ 

with  ei  =  (1, 0, . . . ,  0). 


2-? 


6n{x1)'Z2{x2)  •  •  .'f' d{xd)dx , 


Let  K0  =  5(0,1)  and  if  /  >  0,  Kt  =  5(0,  2l)  -  5(0,  2'”1).  If  *  e  Kh  — V-  and 


x  +  k  ne  i 

21  2.21 


2-? 


belong  to  B(x o,  |*o  —  A«a|  +  n2l  1).  Therefore 

Jk  (Alei/2/)  6>n(aq)^2(a;2)  .  ■  ■  U(xd)Gfa 

<cf  (ko  ~  |  +  n2l-J)a  ^  ^  (|a?o  ~  Pa|  +  n2l~^)c 


Ik,  (1  +  |*|)“+d+1 


2;(«+i) 


which  is  bounded  by  C 


|*o  -  A* a| 


2-oj 


2;(a+i)  if  n2(  J  ^  1*0  -  I,  and  else  by  .  The 


sum  over  l  is  therefore  clearly  bounded  by  |a?o  —  A*a|“  +  2_"J;  hence  (2.14)  holds 
and  the  converse  part  in  Proposition  6  is  now  a  consequence  of  the  converse  part 
of  Theorem  1. 


The  following  corollary  follows  from  Theorem  1  and  Proposition  6.  It  character¬ 
izes  the  Holder  exponent  by  local  decay  conditions  of  the  oscillation,  or,  equivalently, 
by  local  decay  conditions  of  local  suprema  of  the  wavelet  coefficients. 


Corollary  1.  If  /  is  uniformly  Holder,  the  Holder  exponent  of  /  can  be 
computed  using  either 


(2.23) 
or 

(2.24) 


hf(x  o)  =  liminf 

j-t  °o 


\og{dj(x0))\ 
log(2-l)  J  ’ 


hf(x  o)  =  sup 
nets 


lim  inf 

e->0 


l°g(|QS/ (5(a;0,  e))|) 
loge 


Remark:  The  supremum  over  n  in  formula  (2.24)  should  be  understood  as 
follows:  The  term  between  brackets  is  equal  to  n  as  long  as  n  <  [hf(x o)]  and 
becomes  independent  of  n  if  n  >  [hf(x o)]. 
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2.6.  Notes.  The  first  construction  of  an  orthonormal  wavelet  basis  was  per¬ 
formed  by  Haar  in  1909;  we  saw  that  the  corresponding  wavelet  is  discontinu¬ 
ous.  The  first  arbitrarily  smooth  orthonormal  wavelet  bases  were  constructed 
by  Stromberg  in  1981.  Wavelets  in  the  Schwartz  class  were  introduced  by  P.-G. 
Lemarie  and  Y.  Meyer  in  1986.  The  construction  by  multiresolution  analysis  was 
performed  by  S.  Mallat  and  Y.  Meyer  in  1989,  and  compactly  supported  wavelets 
were  introduced  by  I.  Daubechies  shortly  afterwards. 

The  wavelet  characterization  of  the  space  C"(Md)  is  remarkable;  it  is  in  sharp 
contrast  with  Fourier  expansions:  Indeed,  the  space  of  periodic  Ca  functions  cannot 
be  characterized  by  a  condition  bearing  on  the  moduli  of  their  Fourier  coefficients, 
see  [102], 

Hardy  and  Littlewood  are  the  forerunners  of  the  wavelet  method  used  to  prove 
irregularity  at  a  point.  They  had  noticed  that,  if  /  is  smooth  at  xo,  the  convolution 
of  /  with  a  smooth,  well  localized  function  ip  of  vanishing  integral  (a  “wavelet”)  has 
to  be  small  near  xo,  see  [40,  41].  Irregularity  results  are  obtained  by  contraposition. 
(One  can  interpret  wavelet  coefficients  as  a  convolution  product  of  f(x)  with  ip(2J  x), 
this  convolution  being  then  sampled  at  the  points  k 2~J .)  Hardy  and  Littlewood 
applied  this  method  to  the  Weierstrass  functions,  but  also  to  the  “non-differentiable 
Riemann  function”  JP  sm(2™  ^ ,  which  is  much  more  difficult  to  handle  since  its 
Holder  exponent  is  everywhere  discontinuous;  they  showed  in  1914  that  this  function 
is  not  differentiable,  except  perhaps  at  the  rational  points  ^q+i  >  and  J-  Gerver 
showed  in  1970  that,  indeed,  it  is  differentiable  at  these  points,  see  [37].  This  was 
the  first  function,  with  unbounded  variation,  which  was  proved  to  be  multifractal 
(in  1996,  see  [51]).  More  properties  of  this  remarkable  function  are  exposed  in  [62]. 

If  irregularity  results  using  the  wavelet  technique  were  at  least  implicit  in  Hardy 
and  Littlewood’s  paper  in  the  1910s,  the  first  converse  result  only  appeared  in 
1989,  in  [47];  nonetheless,  one  should  mention  a  remarkable  result  of  J.-P.  Kahane 
who  proved  in  1976  the  existence  of  slow  points  of  the  Brownian  motion  using  its 
decomposition  on  the  Schauder  basis,  which  predates  these  wavelet  techniques,  see 
[63,  64], 

We  mentioned  in  Section  1.6  several  generalizations  of  pointwise  Holder  regular¬ 
ity.  Another  generalization,  the  scaling  exponent,  is  based  on  the  Littlewood-Paley 
decomposition  of  /;  it  has  been  introduced  by  Y.  Meyer  in  [88].  Moduli  of  conti¬ 
nuity  more  general  than  6(x)  =  xa  and  their  wavelet  characterization  are  studied 
in  [60]. 

The  idea  of  studying  the  properties  of  a  function  using  the  wavelet  leaders  d\ 
rather  than  the  wavelet  coefficients  first  appeared  in  the  resolution  of  the  following 
problem:  Find  a  formula  which  yields  the  upper  box  dimension  of  the  graph  of  a 
function;  the  answer  is  expressed  in  terms  of  the  d\,  see  [53]. 

Condition  (2.13)  used  to  be  written  in  the  form 

|cA|  <  C2~aj{l+\2jx-k\)a 

(which  is  clearly  equivalent  to  (2.13)).  This  condition  is  a  particular  case  of  the 
two-microlocal  conditions  C's,s  ( xo )  defined  by 

(2.25)  |cA|  <  C2~sj(l  +  \2jx  -  k\)~s'. 

Two-microlocal  conditions  were  introduced  by  J.-M.  Bony  as  a  tool  in  the  study  of 
the  propagation  of  singularities  of  solutions  of  PDEs.  The  set  of  couples  (s,  s')  such 


WAVELET  TECHNIQUES  IN  MULTIFRACTAL  ANALYSIS 


23 


that  /  G  Cs,s  (xo)  is  called  the  two-microlocal  domain  of  /  at  xq.  It  yields  a  very 
accurate  information  on  the  local  oscillations  of  /  near  xo  and  the  regularity  of  its 
fractional  derivatives  and  primitives  at  xo-  The  properties  of  the  two-microlocal 
domain  were  investigated  by  J.  Levy-Vehel,  S.  Seuret  and  their  collaborators,  see 
[78],  and  the  paper  in  the  present  volume  [79]. 

The  converse  part  in  Theorem  1  requires  a  uniformly  Holder  assumption.  Even 
if  this  assumption  can  actually  be  slightly  weakened,  a  uniform  regularity  assump¬ 
tion,  strictly  stronger  than  continuity,  is  necessary,  as  shown  by  Yves  Meyer,  see 
[48].  Therefore,  wavelet  methods  cannot  be  used  to  study  functions  which  have  a 
dense  set  of  discontinuities.  Such  functions  are  not  just  a  mathematical  curiosity; 
for  instance,  most  Levy  processe  (processes  with  independent,  stationary  incre¬ 
ments),  which  play  a  central  role  in  probability,  share  this  property,  see  [54].  It  is 
also  the  case  for  Davenport  series  which  are  of  the  form 

OO  ^ 

an{nx}  where  {x}  =  x  —  [x]  —  — , 

n  =  1 

see  [58].  Note  however  that,  in  [61],  pointwise  Holder  estimates  are  obtained  with 
only  a  Besov  regularity  assumption  which  implies  no  uniform  Ce  regularity  (but 
is  strictly  stronger  than  continuity).  Similarly,  the  logarithmic  correction  in  the 
converse  part  of  Theorem  1  is  best  possible,  see  [48]. 

The  generalizations  of  the  pointwise  regularity  conditions  mentioned  in  Sec¬ 
tion  1.6  also  have  wavelet  characterizations;  for  instance,  C.  Melot  has  obtained  a 
wavelet  characterization  of  the  T^(x o)  condition,  see  [85]. 

The  first  quasi-sure  result  concerning  Holder  regularity  is  a  famous  theorem  by 
Banach,  in  1931,  which  states  that  a  continuous  function  is  almost  surely  nowhere 
differentiable,  see  [13].  Z.  Buczolich  and  J.  Nagy  proved  in  [22]  that  quasi-all 
monotone  continuous  functions  are  multifractal  and  their  spectrum  is 

d(H)  =  H  if  H  G  [0,  1], 

=  — oo  if  H  >  1. 

3.  The  multifractal  formalism 

The  spectrum  of  singularities  of  many  mathematical  functions  can  be  deter¬ 
mined  directly  from  its  definition.  On  the  opposite,  for  many  real-life  signals, 
whose  Holder  exponent  is  expected  to  be  everywhere  discontinuous,  the  numerical 
determination  of  their  Holder  regularity  is  not  feasible,  and  therefore,  one  cannot 
expect  to  have  direct  access  to  their  spectrum  of  singularities.  In  such  cases,  one 
has  to  find  an  indirect  way  to  compute  d(H );  the  multifractal  formalism  is  a  for¬ 
mula  which  is  expected  to  yield  the  spectrum  of  singularities  of  /  from  “global” 
quantities  which  are  numerically  computable.  Mathematically,  these  quantities  are 
interpreted  as  indicating  that  /  belongs  to  a  certain  subset  of  a  family  of  function 
spaces.  There  exist  slight  variants  on  these  formulas,  which  depend  on  the  family  of 
function  spaces  used.  Their  validity  never  holds  in  complete  generality.  However, 
three  types  of  verification  can  be  performed: 

•  The  multifractal  formalism  is  proved  under  additional  assumptions  on  / 
(a  self-similarity  assumption,  or  for  a  class  of  particular  random  processes, 
or  even  for  very  specific  functions  /,  see  [3,  50,  51]  for  instance). 
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•  It  is  proved  for  a  “large”  subset  of  the  function  space  considered;  the 
quasi-sure  result  of  Proposition  5  in  Section  2.4  supplies  such  an  example. 
We  will  also  see  prevalence  results  in  Section  5. 

•  The  multifractal  formalism  is  shown  to  yield  an  upper  bound  of  the  spec¬ 
trum  of  singularities  of  any  (uniformly  Holder)  function,  see  Section  3.2. 

In  this  section,  we  will  first  describe  the  heuristic  arguments  used  in  the  deriva¬ 
tion  of  the  multifractal  formalism.  They  will  be  presented  in  two  different  contexts: 
First,  in  Section  3.1  using  wavelet  expansion,  and  in  Section  3.3  directly  on  incre¬ 
ments  of  the  function.  This  will  lead  us  naturally  to  introduce  some  new  function 
spaces,  the  oscillation  spaces,  which  will  be  studied  afterwards.  In  Section  3.2,  we 
prove  that,  indeed,  the  multifractal  formalism  yields  a  general  upper  bound  of  the 
spectrum. 


3.1.  Derivation  of  the  multifractal  formalism:  The  wavelet  method. 

A  multifractal  formalism  is  a  formula  which  is  expected  to  yield  the  spectrum 
of  singularities  of  a  function  from  the  estimation  of  ‘global’  quantities  which  are 
numerically  computable.  The  first  formula,  based  on  Lp  norms  of  the  increments 
of  the  function,  was  proposed  by  G.  Parisi  and  U.  Frisch  and  in  1985,  see  [92]. 
Soon  afterwards,  A.  Arneodo,  E.  Bacry  and  J.-F.  Muzy  proposed  an  alternative 
formula  based  on  the  continuous  wavelet  transform,  see  [4,  7].  In  this  section,  we 
derive  a  slightly  different  formula  based  on  the  wavelet  leaders  d\.  The  advantage 
will  be  that  one  does  not  have  to  assume  that  the  Holder  singularities  are  “cusp- 
singularities” ;  a  cusp  singularity  at  xo  displays  a  behavior  similar  to  \x  —  xo\a,  as 
opposed  to  chirps  which  display  strong  oscillations  like  \x  —  *o|“sin(|*  —  *o|_/3), 
see  [6]. 

Using  heuristic  arguments,  we  now  derive  this  multifractal  formalism.  We  will 
form  ‘global  quantities’  based  on  the  wavelet  leaders  d\.  Indeed,  if  /  is  uniformly 
Holder,  which  we  assume  from  now  on,  Theorem  1  shows  that,  if  hj  ( xo )  =  H ,  there 
exists  an  infinite  number  of  dyadic  cubes  A  which  contain  xo  (or  are  adjacent  to 
such  a  cube)  such  that 


lim 

j-y+o 


loSdx  ) 

l°g(2 -d)) 


which  we  will  write  d\  ~  2~Hi .  Global  quantities  which  are  natural  to  consider 
are  the  /^-averages  of  the  d\  at  all  scales.  In  order  to  keep  as  much  information 
as  possible  on  the  d\,  we  don’t  restrict  the  computation  to  p  >  1  or  even  p  >  0; 
therefore,  if  S2  is  a  bounded  domain  of  and  p  £  M*,  let 


cO  _  c)-dj 


£ 


AeAj-nn 


where  the  sum  is  taken  on  all  wavelet  leaders  of  scale  j  corresponding  to  dyadic 
cubes  localized  in  fh  The  order  of  magnitude  of  S^j  when  j  — >■  +oo  is  estimated 
with  the  help  of  the  scaling  function  of  /  on  S2  defined  by 


( lQg  (SP,j)  ) 

{  l°g(2--0  J  ’ 


(3.1) 


uif  (p)  =  lim  inf 

J  a  — 4.  -1-  rv~ i 
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which  means  that  •  is  of  the  order  of  magnitude  of  2  “/ when  j  —y  +oo. 
More  precisely,  by  definition  of  w  j?  (p) , 

VS  >0,30  0  J2  d\ <  C2Sj 

AeAj-nn 


and  there  exist  jn  —y  +oo  such  that 

VS  >  0,  3C  >  0  J2  dl 2(“/(p)-d)j'“  >  C2~Sj> 


AGAJrenn 

These  two  conditions  can  be  rewritten  as  follows: 

(3-2)  t  E  dPx2^?(p)-d+s)i  — >  +oc 

j=o  xeAjnfi 


For  any  J  >  0, 
when  J  — >  +oo 


and 

(3.3) 


E  E  dpx2^w?^-d-^j  <  C. 

j=o  xeAjnfi 


Let  xo  E  and  e  >  0.  If  hj(x o)  =  H,  there  exists  a  sequence  of  dyadic  cubes  A  in 
the  influence  cone  at  xq  such  that 


(3.4)  |A|  <  e  and  <  dx  <  . 


For  each  xo  E  £2  at  which  the  Flolder  exponent  takes  the  value  H ,  let  us  pick  such  a 
cube  A,  and  denote  by  A (H)  the  set  of  these  cubes  and  their  adjacent  neighbours. 
The  cubes  which  belong  to  A (H)  therefore  constitue  an  e-covering  of  Eh-  For  this 
particular  subset  of  dyadic  cubes,  (3.3)  implies  that 

E  <  c, 

XeA(H) 

and  therefore,  taking  (3.4)  into  account, 

VS  >  0,  E  |A| (-“?(p)+d+HP+Sp+s)  <  c- 

XeA(H) 


thus  d^(H)  (which  denotes  the  spectrum  of  singularities  of  /  restricted  to  £2)  sat¬ 
isfies 

d^{H)  <  uif{p)  -d-  Hp. 

Furthermore,  (3.2)  implies  that,  for  at  least  one  value  of  H ,  which  we  will  denote 
by  H0, 

E  dpx2{^{p'>-d'>j  >  C2~Sj 

XeA(H0) 

(because  otherwise,  when  one  ‘sums  over  all  values  of  H’,  the  corresponding  quan¬ 
tity  would  decrease  exponentially).  Thus  there  exists  H o  such  that 

d^f(Ho)  =  —  (p)  +  d  +  Hop- 

Summing  up,  we  see  that  the  spectrum  of  /  restricted  to  £2  satifies,  for  any  p  0, 
d*f{H)  <  -uif(p)  +  d  +  Hp  and  3H0  :  d^(H0)  = -uif(p)  +  d  +  H0p; 

which  can  be  rewritten  as  follows: 


(3.5) 


u/}{p)  =  inf(d  —  d%(H)  +  pH). 

J  H  J 
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In  order  to  determine  the  spectrum  on  all  of  Md,  writing  as  a  countable  union 
of  bounded  domains,  we  have 

df(H)  =  sup 
ncKd 

and  therefore,  we  have  to  replace  in  (3.5)  the  function  ui^(p)  by 

(3.6)  cjf(p)  =Mcjf(p). 

We  thus  obtain 

(3.7)  cjj(p)  =  mi(d-  df{H)  +pH). 

ti 

We  will  show  in  Section  4.3  that  ojf(p)  (which  is  defined  on  M\  {0})  extends  to 
a  concave  function  on  M.  Denote  by  df(H)  the  concave  hull  of  df(H)  (i.e.  the 
smallest  concave  function  larger  than  df(H)).  Formula  (3.7)  can  be  interpreted  as 
stating  that  —df(H)  and  —uif(p)  +  d  are  convex  conjugate  functions,  and  each  can 
therefore  be  deduced  from  the  other  by  a  Legendre  transform.  It  follows  that 

(3.8)  df(H)  =  inf  (Hp  —  ujt(p)  +  d)  . 

p 

Formula  (3.8)  tells  us  that  we  only  expect  to  recover  the  concave  hull  of  the  spec¬ 
trum  of  singularities  from  the  scaling  function.  However,  it  often  happens  that 
df(H)  actually  is  a  concave  function,  in  which  case  (3.8)  allows  us  to  recover  df(H) 
completely.  When  such  is  the  case, 

(3.9)  df(H)  =  inf  (Hp  -  uif(p)  +  d)  , 

p 

and  we  will  say  that  the  multifractal  formalism  is  satisfied  by  /. 

3.2.  Upper  bound  of  the  spectrum  of  singularities.  When  p  is  positive, 
the  scaling  function  can  be  given  a  functional  interpretation.  For  that,  we  introduce 
the  oscillation  spaces  which  are  defined  by  the  following  condition  on  the  wavelet 
expansion. 

Definition  15.  Let  s£l  and  p  >  0;  a  distribution  /  belongs  to  CF(Md)  if  its 
wavelet  coefficients  satisfy 

J2n\p<c 

k 

and 

(3.10)  30  0  Vj  >  0  2^-^  J2  <  C. 

k 

Remarks:  These  function  spaces  are  a  particular  case  of  the  spaces  Osp’s  (Md) 
previously  considered  in  [53,  57]. 

It  follows  immediately  from  Definition  15  that,  if  p  >  0,  then 

0Jf{p)  =  sup{s  :  /  £  06pffoc}, 

therefore  the  function  ojf(p)  indicates  which  Op  spaces  contain  the  function  /. 

One  easily  checks  that  the  oscillation  spaces  are  Banach  spaces  if  p  >  1;  they 
are  complete  metric  spaces  if  0  <  p  <  1  for  the  metric  defined  as  follows:  Let  / 
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and  g  be  two  functions  in  Ow .  The  coefficients  of  /  are  denoted  c\  and  Ck  as  in 
Definition  15,  and  those  of  g  are  denoted  c\  and  Ck]  then 


(3.11) 

d(f,g)  =||  f-g  || 


W  | Ck  -  ck\p  +  sup 

V  j>° 


(  2{sP~d'>j  V  (  sup  |ca< 

V  k  VvcA 


c\> 


The  following  result  states  that  the  multifractal  formalism  yields  an  upper 
bound  for  the  spectrum  of  any  function. 


Theorem  2.  Let  /  be  a  uniformly  Holder  function.  Then 
(3.12)  df{H)  <  inf  (Hp-uf{p)  +  d)  . 

p£M  * 

First,  we  consider  the  case  where  p  is  positive.  Recall  that  the  sets  Gh  were 
defined  in  (1.5). 


Proposition  7.  Let  p  >  0;  if  /  is  a  uniformly  Holder  function  in  Op, 

(3.13)  V.H  >  s - ,  dim^Gn)  <  d  +  Hp  —  sp. 

Furthermore,  if  s  —  |  >  0,  then  Gh  =  0  for  any  H  <  s  — 


Proof:  With  regards  to  the  case  H  <  s  —  if  /  E 

|cA|  <  dx  <  C2-(s-^j 

so  that  V*o,  hf(xo)  >  s  —  j. 

Let  us  now  prove  the  first  assertion.  Let 

GjtH={  A:  Ma|  >  2-^}  , 

and  denote  by  Njtn  the  cardinality  of  Gjtn-  By  hypothesis,  /  £  O'  so  that 
2(°P-d)iJ2\d\\p  <  C,  and 

2(‘P-d)jNjiH2-Hpj  <  c, 
so  that  NjtH  <  C2^-sp+hp^. 

Denote  by  Fjh  the  set  of  cubes  A  of  scale  j  such  that  either  A  £  Gjtn,  or  A  is 
adjacent  to  a  cube  of  Gjtn-  Clearly, 

Card(Fj  H)  <  3 dCard{GjiH)  <  3dC2 (d-sp+HP)j . 


Denote  by  Fh  =  limsupFj  iy  the  set  of  points  that  belong  to  an  infinite  number  of 

i->+oo 

Fj,H-  If  xo  (fi  Fh,  then  there  exists  jo  (=  jo(x))  such  that,  for  any  j  >  jo,  we  have 
dj(x o)  <  2~h^\  thus  we  can  choose  C  (  =  C(x))  large  enough  so  that 

Mj  >  0,  dj(x0)  <C2~Hj . 

Theorem  1  implies  that  /  £  Cfl  (x o);  so  that  Gh  C  Fh- 

It  remains  to  bound  the  dimension  of  Fh-  Let  e  >  0,  and 

jo  =  inf{j  :  \[d2~3  <  e). 
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We  choose  for  e-covering  of  Fh  all  the  cubes  A  such  that  j  >  jo  and  A  E  Fj,H- 
Clearly, 

00  s 

Y  Diam(B\)s  <  C  Y  Card(Fj)H)  (Vd 2“^ 

j=3a 

2(d-sp+Hp-6)j 

3=3  0 

which  is  hnite  if  i  >  d  +  Hp  —  sp ;  hence  the  first  part  of  the  proposition  holds. 

Let  us  now  check  that  the  last  point  of  the  proposition  follows.  If  xo  E  Eh, 
hf(x 0)  =  H,  and  Vi?'  >  H,  x 0  E  G#*;  so  that  Eh  C  G#o  Let  p  >  0;  by  definition 
of  ujf(p)  we  have  Ve  >  0,  /  E  ojiF  e^P ,  so  that 

MH '  >  Ed ,  Me  >  0  d(-ff)  =  dimiEn)  <  dim(GH')  <  rf  +  E['p  —  uj j(p)  +  e, 
and  thus  the  spectrum  of  /  restricted  on  any  bounded  set  SI  satisfies 

dn(H)  <  d  +  Edp  —  uif(p). 

Since  this  upper  bound  is  valid  for  all  p  >  0,  (3.12)  follows  (with  the  inhmum  taken 
only  on  M+). 

We  consider  now  the  case  where  p  is  negative.  In  this  case,  we  will  obtain  a 
result  which  is  stronger  than  Theorem  2  since  it  yields  a  bound  for  the  packing 
dimension  of  the  Holder  singularities. 

Proposition  8.  Let  p  <  0;  if  /  E  Osp,  the  packing  dimension  of  Bh  is  bounded 
by  d  —  sp  +  Hp. 

Note  that,  in  contradistinction  with  the  case  p  >  0,  we  do  not  have  to  make 
any  uniform  regularity  assumption  here. 

Proof:  Let  S  >  0  and  J  such  that  2-J  <6.  If  /  E  C'H(x 0),  there  exists  A  >  0 
such  that 

Vj  >  J,  sup  |cA'|  <  A2~Hj-, 

A'C  3,j(xo) 

so  that,  since  p  <  0, 

(3.14)  f  sup  \cx,\)  >  Ap2~HpE 

yA'CAjfrco)  J 

Denote  by  f l a  the  set  of  points  x  where  (3.14)  holds  for  any  j  >  J .  Clearly, 

Bh  C  IU, 

T>0 

where  the  union  can  be  written  as  a  countable  union.  Since  /  Ed),  there  are  at 
most  CAP 2A~sp+Hp+f)3  cubes  A  satisfying  (3.14),  so  that  the  upper  box  dimension 
of  each  set  f l a  is  bounded  by  d  —  sp  +  Hp.  The  proposition  follows  by  countable 
union,  and  the  upper  bound  (3.12)  also  follows  (with  the  inhmum  taken  on  M-); 
therefore  (3.12)  is  completely  proved. 
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3.3.  Derivation  of  the  multifractal  formalism:  The  oscillation  method. 

In  order  to  derive  a  multifractal  formalism,  we  started  from  (2.23);  however,  one 
can  alternatively  start  from  the  characterization  of  pointwise  regularity  based  on 
local  oscillations  (2.24):  If  hf(x o)  =  H,  there  exists  an  infinite  number  of  balls 
B(x o,  £j)  such  that,  for  n  large  enough, 


lim 

c  j  — y  0 


log  (|OS” c y ) ) | ) 

log(ei) 


=  H. 


In  this  formulation,  the  “global”  quantities  which  are  natural  to  consider  are 


(3.15)  T^e  =  edj  \OSnf(B(l,3Vde))\p 

l£el,d  Of! 


(we  use  balls  of  radius  3 \fde  so  that  there  is  enough  overlapping,  but  the  precise 
value  3  is  not  important).  Let 


(p)  =  lim  inf 


d+lim  inf 

e-MI 


^logf  X)  \OSnf(B(l,3y/de))\^^ 

yeezdnn 

log(e) 


The  argument  we  developed  using  the  can  be  reproduced.  In  the  present 
setting,  it  leads  to  the  following  formulation  of  the  multifractal  formalism: 


(3.16)  d(H)  =  inf  (Hp  —  vt(p)  +  d) , 

p 

where 

VJ{P)  = 

When  p  is  positive,  the  function  vj{p)  also  has  a  function  space  interpretation. 

Definition  16.  Let  p  >  0  and  s  £  M.  A  function  /  belongs  to  Vp(Md)  if 
/  E  L°°(lStd)  and  if,  for  n  >  [s], 

(3.17)  300  Ve  >  0  ed~sp  \OSnf(B(l,  3vQe))\p  <  C. 

Note  that 

(3.18)  u/(p)=sup{S:/eV;/p0C(Md)}. 

The  function  vj  (p)  indicates  to  which  V‘  spaces  the  function  /  belongs  locally. 

Before  considering  the  problem  of  the  validity  of  the  multifractal  formalism, 
one  should  first  check  that  its  two  formulations  using  either  wavelets  or  oscillations 
lead  to  the  same  formula.  There  exists  only  a  partial  result  in  this  direction  which 
concerns  positive  values  of  p.  We  will  show  that,  if  /  is  uniformly  Holder,  for  any 
P>  0, 

^f(p)  =  Vf(p)- 

This  result  is  a  consequence  of  the  following  theorem. 
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Theorem  3.  The  following  imbeddings  hold 
(3.19)  Ve  >  0,  V?7  >  0  Ce  n  Vs+V’p  ^  Osp  ^  Vs’p; 

thus,  if  /  :  — >■  M  is  uniformly  Holder,  then  the  p-oscillation  exponent  of  /  is 

given  by 

uif(p)  =  sup{s  :  /  E  Op/p’°} 


(3.20) 


log  X) 


lim  sup  ■ 

i->+oo 


,  AGA, 


sup  |CA' 

A'CA 


j  log  2 


Let  us  stress  the  fact  that  this  result  deals  only  with  the  case  p  >  0.  We  do  not 
know  if  the  scaling  function  u>f(p)  can  be  characterized  by  local  oscillations  when 
p  is  negative. 


3.4.  Proof  of  Theorem  3.  We  assume  that  we  use  compactly  supported 
wavelets.  (The  general  case  will  follow  afterwards  using  Theorem  4  which  states 
that  the  scaling  function  is  independent  of  the  wavelet  basis).  Coming  back  to  the 
comparison  of  wavelet  coefficients  and  oscillations  performed  in  Section  2.5,  one 
sees  that  the  function  9n  (a?i ) SP2 (a?2 )  •  •  .  xf,c!(*d)  which  appears  in  (2.22)  is  compactly 
supported,  say  in  [—  A,A\d.  Let  B\  =  B(k2~i ,  A\fd2~i).  It  follows  from  (2.22)  that 
|cA|  <  C  •  OS "(5a);  thus,  if  A'  C  A, 

\cx,\  <  C  ■  OS]  (Bx,)  <  C  ■  OS]  (Bx); 
so  that  d\  <  C  ■  OS "(5A).  Therefore, 

2(sp-d)j  J2  d{  <  C2<-sp~d',j  J2  [°s/(5a)]P  • 

A  G  A  j  A  G  A  j 


In  order  to  obtain  a  converse  estimate,  we  first  assume  that  /  is  compactly 
supported.  We  will  bound  the  oscillations  of  each  component  fji  on  the  ball  B\  by 

A£  A  j 

If  j1  >  j,  then 

fjl  =  ^2  cx’^  A'- 

A'GA,/ 

Using  the  localization  of  the  wavelets,  it  follows  that 


OSnf.,(Bx)  <C\\f\\L~{Bx)<C  sup  |cA-|. 

3  CA'nBA^0 


In  particular,  since  /  E  C  e(Md),  it  follows  that 

(3.21)  OS’}  (Bx)  <  C2~ej' . 


Let  d\  =  sup  sup  \Cp,\.  We  have 
C  An^^0  /ii'C/ii 


(3.22) 


OS  ll(Bx)<Cdx. 
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Since  d\  is  a  supremum  of  d\  on  at  most  (C  +  2)d  cubes,  it  follows  that 
(3.23)  £  A<(C  +  2U  A- 

A  £  A  j  A  G  A  j 

If  j'  <  j,  applying  Taylor’s  formula  to  A ,  it  follows  that,  if  x  and  x  +  nh  belong 
to  B\, 

\Ahfj'  I  <  c\h\n  sup  \dafj'(x)\, 

xeB  A 

where  the  supremum  bears  on  the  |a|  =  n.  But,  if  |a|  =  n, 

\dafj-{x)\<C2nj  sup  |cA'|. 

CA'nBA7£0 


Since  \h\  <  C 2  J,  it  follows  that 

OSr,(Bx)  <C2-Ci-i')  sup  |CA,|. 

3  C\'nBx7Bt) 

Therefore,  if  AJ  denotes  the  dyadic  cube  of  scale  j '  which  includes  A, 

(3.24)  OS’}.,  (B\)  <  C2-n('j-j'^dX]l . 

OO  OO 

Since  /  =  fji,  we  have  OS "(_BA)  <  OS".;(5A).  Using  (3.21),  it  follows 

i'=-i  i'=-i 

that 


and  using  (3. 22), we  get 


£  OS^,(5A)<U2-^2 


£0S?  (BA)<Cj2dA. 


Therefore 


Since 


it  follows  that 


j- 1 

OS ](BX)  <  £  OS/y  (5a)  +  Cj2dx  +  C' 2- 
i'=-i 


Vp  >  0,  I  £  a,  j  <  (n  sup  a,)p  <  np  £(ai)? 


vi=i 


8  =  1 


i-i 


(OS" (B\))p  <  C(j  +  3)  £  [0Snfi,(Bx)]P +j2P(dx)P +  2~^2 


\j'=- 1 

Thus,  using  (3.24),  it  follows  that 

j- 1 


(OS" (Bx))p  <  C(j  +  3)  £  2 -"PU-j')  (dxj,y  +  fP(dx)P  +  2 


-epj 


U'=-i 
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When  we  sum  on  A,  each  term 


Xj/  j  appears 


2dU  3  )  times,  so  that 


E  (°s/(5a))p  < 

AG  A  j 


cu  +  3)(  E  2(d~np)ii~i')  E  (v)P +fp  E  (Aop  +  C2dj2~epj2 ) 

y'=-i  A'GAj;  ag  Aj  / 

(in  the  derivation  of  the  last  term  we  use  the  fact  that  /  is  compactly  supported, 
so  that  the  sum  bears  on  at  most  C 2dj  dyadic  cubes  of  scale  j).  Using  (3.23),  we 
can  replace  d\  by  d\.  By  hypothesis, 

V/  Y,  (rfA'f  <  C2(d-‘P+£V'-, 

A'GAj; 

so  that 

E  (os/(5a))p 

AgA  j 


(  \ 

<  c(j  +  3)  E  2-(np-d)(j-j')2(d-sp  +  e)j'  _|_  j2p  2<d- sp  +  e)j  _|_  (^7  2  2  —  2  1 

v--1  / 

which  is  bounded  by  2C-sP+2eb  if  n  is  picked  large  enough  so  that  np  —  d  >  sp  —  d, 
i.e.  n  >  s. 

The  result  follows  for  any  function  /  (non  necessarily  compactly  supported), 
using  a  partition  of  unity:  Let  9  be  a  C°°  compactly  supported  function,  satisfying 
9(x  —  k)  =  1.  We  apply  the  result  to  each  fk(x)  =  f(x)9(x  —  k)  and  use 

k  GZ 


II  3  I  E  II  Ac  II  o°p  • 

k 

A  similar  argument  also  works  for  the  V* . 

Let  us  now  check  that  the  £ios;  (B(p\,  AA))  |P  are  independent  of  the  (large 
enough)  constant  A.  We  will  check  that,  if 

A=  E  los"(3A)f 

AG  A  j 


and 

Bj  =  E  |OS^(AA) 

AgA  j 


where  A  >  3,  then 


Of  course,  Aj  < 


=  liminf 
i->  +  oo  log(2  J)  j-¥  +  OO 

Bj.  Let  l  be  such  that 


lQg  (Bj) 
log(2-J)  ’ 


^2l  <  A  <  2l . 
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Then 

Bj  <  Y  |OS^(2'A)|P  <  3d2dl  Y  |OS?(3A')|P 

AgA  j  X1  £  A  j-i 

because  each  cube  of  scale  j  is  contained  in  a  cube  3A',  where  A'  is  a  dyadic  cube 
of  scale  l  —  j,  and  each  cube  of  this  type  contains  at  most  3d2dl  cubes  of  scale  j . 

In  order  to  end  the  proof  of  Theorem  3,  we  have  to  check  that 

ed~‘p  Y  IOS?  {B{l,3Vdc))\p  <  C2^p-d^j  Y  los/UA)r 

JG  AG  A  j 

wher  we  can  pick  any  A  >  3.  we  pick  j  such  that 

2~j  <  e  <  2  •  2_J , 

and  A  large  enough  so  that  any  ball  B((l,  3 y/de)  is  entirely  included  in  a  cube  AX 
(A  G  A  j).  The  estimate  follows. 

3.5.  Notes.  We  mentioned  that  the  multifractal  formalism  may  not  yield  the 
correct  result  when  applied  to  chirps;  this  problem  is  discussed  in  [6]  and  [86].  A 
chirp  at  xo  displays  strong  oscillations  near  xo  similar  to  those  of 

\x  —  x0\a  sin 

Their  study  was  initiated  by  P.  Tchamitchian  and  B.  Torresani  in  [98],  see  also 
[24];  there  exist  two  slightly  different  ways  to  model  such  such  a  behavior,  either 
mathematical  chirps,  see  [60,  88],  or  oscillating  singularities,  see  [6,  8].  In  both 
cases,  the  couple  (a,  f3)  is  used  to  label  either  chirps  or  oscillating  singularities. 
Therefore,  one  can  associate  to  such  singularities  a  new  spectrum  d(a,  /?)  indexed  by 
two  parameters  which  is,  as  expected,  the  Hausdorff  dimension  of  these  singularities. 
Such  two-variables  spectra  have  been  determined  either  for  stochastic  processes 
(random  wavelet  series),  see  [9]),  or  for  quasi-all  functions  in  Besov  or  Sobolev 
spaces,  see  [86]. 

The  justification  of  the  multifractal  formalism  that  we  proposed  is  an  adapta¬ 
tion  of  the  one  initially  introduced  by  G.  Parisi  and  U.  Frisch  in  [92]  (see  also  [39]). 
One  difference  is  that  the  “local  quantities”  which  were  their  starting  point  differ 
from  the  Tp  e:  Let 

SP{  0=  J \f{x  +  l)~f{x)\Pdx; 

the  corresponding  scaling  function  used  in  the  former  formalism  is 
C/ (p)  =  sup{r  :  \Sp(l)\  <  C\l\T}- 

Recall  that,  if 

(3.25)  dQ-1^  <  s  <  1, 

the  Besov  space  Bsp’q  can  be  characterized  by  the  condition 

[  sup  (  f \f(t  +  h)  -f(t)\pdt)  -^-<G, 

JO  0<|/i|<t  \J  J  1 

see  [100],  This  characterization  clearly  implies  that  the  corresponding  scaling  func¬ 
tion  is 

(3.26) 


C/  (p)  =  sup{s  :  feB’J™}, 
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at  least  as  long  as  d  ^  <  1  (we  can  choose  any  q  >  0);  in  general, 

one  should  rather  take  (3.26)  as  a  definition  of  the  scaling  function.  One  easily 
checks  that  Q  is  also  concave.  The  corresponding  multifractal  formalism  asserts 
that 

(3.27)  df(H)  =  inf  ( Hp  -  Q  (p)  +  d)  . 

p 

We  will  show  in  the  following  section  (Proposition  9)  that,  if  s  >  j,  Op  =  Bp00 . 
It  follows  that,  when  the  inhmum  in  (3.27)  is  attained  for  values  of  p  such  that 
C(p)  >  d,  this  formula  coincides  with  (3.8).  However,  the  formulation  of  the  mul¬ 
tifractal  formalism  based  on  oscillation  spaces  has  a  wider  range  of  validity,  see 
[6].  In  the  applications,  one  prefers  to  use  wavelet  coefficients  in  the  formulation 
of  the  multifractal  formalism,  rather  than  increments,  which  are  numerically  less 
stable;  indeed,  increments  (or  oscillations)  are  very  sensitive  to  noise;  by  contrast, 
since  the  wavelet  coefficients  are  averages  (integrals  against  smooth  functions),  they 
are  numerically  much  less  sensitive  to  noise.  Note  also  that  the  extension  of  the 
corresponding  scaling  function  Cf(p)  to  p  <  0  has  no  theoretical  backing,  whereas 
we  will  see  in  the  next  section  that  the  scaling  function  based  on  oscillation  spaces 
extends  in  a  ‘canonical’  way  to  p  <  0  (see  Section  4.2).  Wavelet  methods  for  the 
numerical  computation  of  spectra  of  singularities  were  developed  by  A.  Arneodo, 
E.  Bacry,  J.-F.  Muzy  and  their  collaborators,  see  [4,  5,  7], 

The  validity  of  the  multifractal  formalism  for  functions  has  been  the  subject  of 
many  papers.  We  won’t  give  a  detailed  review  of  this  topic;  however,  let  us  just 
mention  that  this  validity  has  been  proved  under  assumptions  of  selfsimilarity  for 
the  function,  either  exact  [50],  approximate  [3,  20]  or  statistical  [54].  These  results 
followed  a  similar  line  of  research  for  multifractal  measures.  We  will  see  in  Section 
5  how  generic  results  of  validity  can  be  obtained. 

Assuming  that  the  multifractal  formalism  holds,  a  function  /  will  be  monofractal 
if  and  only  if  its  scaling  function  is  affine.  This  led  to  a  confusion  and  many 
misunderstandings  between  mathematicians  and  signal  processers:  In  the  applied 
community,  a  multifractal  function  is  usually  defined  as  a  function  /  whose  scaling 
function  (either  uif  or  Q)  is  not  an  affine  function.  Needless  to  say,  such  a  definition 
by  the  negative  insures  that  most  signals  met  turn  out  to  be  multifractal...  This 
partly  explains  the  flourishing  literature  on  multifractality  which  has  spread  to  a 
large  numbers  of  applied  fields. 

4.  Properties  of  the  scaling  function 

Our  purpose  in  this  section  is  to  derive  some  properties  of  the  scaling  function 
that  have  been  used  in  the  derivation  of  the  multifractal  formalism.  The  multifractal 
formalism  given  by  (3.9)  relates  the  scaling  function  (which  has  been  constructed 
with  the  help  of  wavelet  coefficients)  to  the  spectrum  of  singularities  (which  has 
been  defined  without  any  reference  to  a  particular  wavelet  basis).  Therefore,  an 
obvious  consistency  requirement  is  to  check  that  the  scaling  function  is  independent 
of  the  choice  of  the  wavelet  basis. 

Theorem  4.  Let  /  £  L°° .  If  the  wavelets  used  belong  to  the  Schwartz  class, 
then,  for  any  p  /  0,  the  scaling  function  of  /  is  independent  of  the  wavelet  basis 
chosen. 


WAVELET  TECHNIQUES  IN  MULTIFRACTAL  ANALYSIS 


35 


The  proof  of  Theorem  4  will  be  split  into  two  cases,  depending  on  whether 
p  is  positive  or  negative.  When  p  is  positive,  we  already  saw  that  the  scaling 
function  indicates  which  oscillation  spaces  /  belongs  to.  We  will  show  in  Section 
4.1  that  their  definition  is  independent  of  the  wavelet  basis  chosen,  thus  proving 
Theorem  4  when  p  is  positive.  The  case  p  <  0  will  be  examined  in  the  Section  4.2. 
The  function  space  interpretation  of  the  scaling  function  is  important  for  another 
reason:  One  cannot  expect  (3.8)  or  (3.9)  to  hold  for  any  function;  indeed,  it  is 
very  easy  to  construct  counterexamples  to  the  multifractal  formalism;  the  central 
question  of  multifractal  analysis  is  to  determine  the  domain  of  validity  of  (3.9).  Up 
to  now,  most  results  concerned  specific  functions  or  random  processes.  A  way  to  see 
how  general  results  of  validity  can  be  obtained  is  to  remark  that,  since  the  scaling 
function  for  p  >  0  implies  that  /  belongs  to  a  certain  intersection  of  oscillation 
spaces,  we  can  wonder  if  there  is  a  “generic”  result  of  multifractality  that  would 
hold  for  ‘most’  functions  of  this  function  space.  We  will  address  this  problem  in 
Section  5. 


4.1.  Oscillation  spaces.  Oscillation  spaces  were  introduced  in  Definition  15; 
they  are  closely  related  to  Besov  spaces.  Let  us  recall  the  definition  of  Besov 
spaces.  If  s  is  large  enough,  Besov  spaces  can  be  defined  by  conditions  on  the  finite 
differences  A f  already  considered  in  Section  2.5:  Let  p,  q  and  s  be  such  that 
0  <  p  <  +oo,  0  <  q  <  +oo  and  s  >  d(^  —  1)+;  then  /  £  Bp’q  if  /  £  IP  and  if,  for 
M  >  s, 

r*  (if 

(4-1)  /  sup  ||  (Af  /)  HI  — T  <  C. 

Jo  0<|ft|<t  1 

The  following  wavelet  characterization  (which  can  be  taken  as  a  definition)  has 
been  proved  in  [87]. 

Definition  17.  Let  s£l  and  p,  q  >  0.  A  distribution  /  belongs  to  B‘’q  if 
the  sequence  cj,  belongs  to  lp  and  if 

•2)  £( 

j> o  \ 


J2  cx2 


<  Too. 


Remark:  This  characterization  is  of  course  reminiscent  of  (3.10);  indeed,  if 
q  =  Too,  /  £  Bp’°°  if  (cfe)  £  lp  and  if 

(4.3)  Vj  ^sp~d^J2\Cx\P  <°o- 

k  ,i 

Since  |ca|  <  |c?a|,  it  follows  that  Op  5*’°°. 

When  p  >  1,  Besov  spaces  are  also  closely  related  with  the  Sobolev  spaces. 


Definition  18.  Let  s  >  0  and  p  >  1.  A  function  /  belongs  to  the  Sobolev 
space  Lp,s(lStd)  if  /  £  Lp  and  if  (Id  —  A)^2  /  £  IP ,  where  the  operator  (Id  —  A)s/2 
is  defined  as  follows:  g  =  (Id  —  A)s/2/  means  that 

v(0  =  (i  +  kl2)s/2/(0 

(the  function  (1  T  l^l2)8/2  being  C°°  with  polynomial  increase,  (1  T  |£|2)s/2/(£)  is 
well  defined  if  /  is  a  tempered  distribution). 
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This  definition  amounts  to  say  that  /  and  its  fractional  derivatives  of  order  at 
most  s  belong  to  Lp ,  see  [87].  If  s  >  0  and  p  >  1,  then 

Bs/  ^  Lp’s  ^  Bsp’°° . 

Thus  Bp,q  is  very  close  to  Lp,s .  The  following  proposition  shows  that  oscillation 
and  Besov  spaces  sometimes  coincide. 


Proposition  9.  If  s  >  |,  then  Op  =  5*’°°. 

Proof  of  Proposition  9:  We  already  saw  that  Op  Bp’°°  .  In  order  to  prove 

the  converse  embedding,  note  that  dpx  =  sup  |ca'|p  can  be  bounded  by  N  N  lcA'|P 

A'CA  CC-  6 

3  >3  A'CA 

(where  the  sum  over  A'  is  taken  on  the  subcubes  of  A  of  scale  j'),  so  that 


J2dPx‘2{sp~d)j 


k  j’>j  A'CA 


sp-d)j 


=  E  icAT2(sp-d)j 

3 '>3  k< 


(4.4)  =  £  £  \cx'\p2{sp~d)j'2~{sp~d){j'~j). 

j‘>3  k' 

If  s  >  y  and  /  G  5p’°°,  then  \c\i  \p2i'sp~d^  <  C;  (4.4)  is  therefore  bounded; 

k‘ 

hence  5*’ 00  Osp. 

The  following  embeddings  will  be  useful  in  order  to  derive  the  properties  of  the 
scaling  function. 

Proposition  10.  (1)  If  si  <  s2,  Osp2  ^  Op1; 

(2)  if  q>P,Osp^OSq~r  +  «; 

(3)  If  /  is  compactly  supported,  then 

(4.5)  /  eOsp  =>  Vg<p,  /  eO). 

Remark:  The  third  embedding  implies  that  Vg  <  p,  Osp  loc  Osq  loc. 

Proof  of  Proposition  10:  The  first  embedding  is  straightforward;  as  regards 
to  the  second  one,  if  /  G  Osp,  for  any  j,  the  sequences  2<'s~p',qd\  are  bounded  in 
lp  uniformly  in  j,  therefore,  they  are  also  bounded  in  lq ,  uniformly  in  j,  if  q  >  p. 
Since 

2  (s-p^d\  =  2  ((s-p+-^--^d\, 

the  second  embedding  is  proved.  As  regards  to  the  third  one,  if  /  is  compactly 
supported,  and  if  the  wavelets  are  also  compactly  supported,  then,  at  each  scale  j, 
at  most  C2dq  of  the  d\  do  not  vanish.  Let  q  <  p;  applying  Holder’s  inequality  to 
the  sequence  dqx,  and  to  the  sequence  identically  equal  to  1  if  d\  0,  and  vanishing 
elsewhere,  we  obtain 

£  d\  <||  d\  ||p/,||  1  ||r, 

AG  A  j 
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where  r  is  the  exponent  conjugate  to  p/q.  This  inequality  implies  that 

q/p 

,  i-a 


Therefore 


ErfA<  Ea]  (c^y 

A  E  A  j  \  A  E  A .j 


J2  2 <  C  |  J2  2 (sP~d)idpx 


q/p 


AE  A  j 


A  E  A  j 


So  that  /  G  Ot  =7  /  E  0\ ’ 


The  spaces  0*  are  dehned  by  conditions  on  the  wavelet  coefficients,  therefore 
the  first  point  to  check  is  that  this  definition  is  independent  of  the  wavelet  ba¬ 
sis  chosen.  In  practice,  one  usually  checks  a  stronger  (but  simpler)  requirement 
which  implies  that  the  condition  considered  has  some  additional  stability;  indeed, 
the  matrix  of  the  operator  which  maps  an  orthonormal  wavelet  basis  onto  another 
orthonormal  wavelet  basis  is  invariant  under  the  action  of  infinite  matrices  which 
belong  to  algebras  M1  that  will  be  dehned  below;  therefore,  one  can  check  that 
Condition  (3.10)  is  also  invariant  under  this  action,  which  is  the  purpose  of  Propo¬ 
sition  11  below;  these  algebras  are  dehned  as  follows: 


Definition  19.  An  infinite  matrix  A(A,  A')  indexed  by  the  dyadic  cubes  be¬ 
longs  to  M1  if 


,  C  2_(^+A)C-f') 

-  (1  +  (j  -  i')2){l  +  2mf(ij')dist(A,  A'))d+7  ' 


Matrices  of  operators  which  map  a  smooth  wavelet  basis  onto  another  one 
belong  to  these  algebras,  and  more  generally  matrices  (on  wavelet  bases)  of  pseu¬ 
dodifferential  operators  of  order  0,  such  as  the  Hilbert  transform  in  dimension  1,  or 
the  Riesz  transforms  in  higher  dimensions,  belong  to  these  algebras,  see  [87].  We 
denote  by  0p(Ad7)  the  space  of  operators  whose  matrix  on  a  wavelet  basis  belongs 
to  M1 . 


Proposition  11.  If  7  >  |s|,  the  operators  which  belong  to  0p(Ad7)  are  con¬ 
tinuous  on  0p(Md). 

Remark:  Since  0p(Ad7)  is  an  algebra,  this  result  clearly  implies  that  the  def¬ 
inition  of  oscillation  spaces  is  independent  of  the  (smooth  enough)  wavelet  basis 
chosen. 


Proof  of  Proposition  11:  Recall  that  we  write  A  =  A  (i,  j,  k),  A'  =  A  (i1 ,  j' ,  k'),... 
see  Section  2.3.  Let  ey  =  A(X' ,  X")c\n  and 


2-(#+7)li- i'l 

WA,A'  “  (1  +  (j  ~  j'PH1  +  2mf(ij')dist(A,  A'))d+7  ' 

The  ey  satisfy  |ev'|  <  C  Cojyyi \cy  |;  assume  that  the  c\  satisfy  (3.10);  we 
want  to  bound 

V 

A"  |CA"  I 
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We  split  the  sum  over  A"  into  two  terms,  depending  on  whether  A"  C  3A  or  A"  3A. 
First  case:  A"  C  3A.  Then, 


z  wA',a"|cA"|< 

A"C3A 


sup  |cA"|  VwA' A"  <  C  sup 

\A"C3A  /  ^77  A"C3A 


CA" 


therefore 


Z  sup  MC 

| A |  =  2  —  J  A'CA  \A"C3A 


P 


<  c  W  sup  |cA//|p 

A"C3A 


<C  3d 

|  A|  =  2-J  +  1 


sup  |cA//|p. 

A"C  A 


Second  case:  A"  <f_  3A.  Let  us  first  prove  that,  in  this  case, 


(4.6) 


wa',a"  <  Cu>xtx"- 


In  order  to  prove  (4.6),  let  us  first  assume  that  j"  <  j.  In  that  case,  dist( A',  A")  > 
^dist( A,  A")  and  \j'  —  j"\  <  | j  —  j"\ ;  (4.6)  follows  directly  from  these  two  results. 

If  j"  >  j ,  then,  since  A'  C  A 


WA'-A"  -  I  (l  +  (/  _  j"Y)dist{\,  A")A+72 lnO' j")(a+7)  ' 
We  now  separate  two  cases: 

If  j  <  j"  <  j',  then 


C  2-(#+7)li-i"l 

iij \ /  \n  <C  - - - 

“  (1  +  (/  -  j"Y)dist( A,  A,,)d+''2f"(d+7)2-(-+A)(i"-i) 
C2-(f+hli-i"l 

“  (1  +  (/  -  i")2)2^f"-fVist(A,  A")d+''2(d+A)i 
C  2 -(#+7)(i"-i) 

-  (1  +  (j  -  j")2)dist{ A,  A")d+'1'2(d+W 

<  Cuix,x"  ■ 


We  finally  consider  the  case  j  <  j'  <  j".  Then 

C  2-(#+7)li-i"l 

- - - 

“  (l  +  (i' -i,,)2)dist(A,A,,)d+''2T(A+7)2-(-+A)(i'-i) 
C  2-(§+7)li-i"l 

“  (1  +  (/  -  i")2)2-^f-f')dist(A,  A")d+''2(d+A)i 
<  Cujx,x"- 

It  follows  from  (4.6)  that 

Z  sup  Z  wa',a"|ca»|  )  <  C  Z  E  WA,A"|cA"| 

|A|  =  2-JA'cA  \X"£3X  J  | A |  =  2  —  A  \  A" 

<  Csup  Z  |cA„|' 

J'  |A'|  =  2-A' 


(the  last  inequality  holds  by  continuity  of  the  operators  in  Op{M1)  on  the  Besov 
spaces  Bp’ 00 ,  see  [87]);  hence  Proposition  11  follows. 
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4.2.  Oscillation  spaces  for  p  <  0.  The  multifractal  formalism  asserts  that 
the  spectrum  of  singularities  of  a  function  /  can  be  deduced  from  its  scaling  function 
u>f(p)  by  a  Legendre  transform.  We  saw  that,  if  p  is  positive,  the  scaling  function 
u>f(p)  determines  which  oscillation  spaces  /  belongs  to.  In  particular,  this  interpre¬ 
tation  shows  that  u>f(p)  is  independent  of  the  particular  wavelet  basis  chosen.  Now 
we  want  to  determine  if  this  interpretation  can  be  extended  to  negative  values  of  p: 
Can  one  define  Op  spaces  for  p  negative?  And,  if  so,  is  the  definition  independent 
of  the  wavelet  basis  chosen?  We  will  need  here  a  different  requirement  than  the  one 
used  in  Proposition  11. 

Definition  20.  An  infinite  matrix  A(A,  A')  is  quasidiagonal  if  A  is  invertible, 
and  if  A  and  A-1  belong  to  MA  for  any  7  >  0. 

The  matrix  of  an  operator  which  maps  a  C°°  orthonormal  wavelet  basis  onto 
another  C°°  orthonormal  wavelet  basis  is  quasidiagonal,  see  [87].  Therefore,  in 
order  to  check  that  a  condition  defined  on  the  wavelet  coefficients  is  independent  of 
the  wavelet  basis  (in  the  Schwartz  class)  used,  one  can  check  the  stronger  property 
that  it  is  invariant  under  the  action  of  quasidiagonal  matrices. 

Definition  21.  Let  C  =  {ca}aga  be  a  collection  of  coefficients  indexed  by  the 
dyadic  cubes.  A  property  V  is  robust  if  the  following  condition  holds:  If  V(C') 
holds  then,  for  any  quasidiagonal  operator  M,  'P(MC)  holds. 

Similarly,  we  will  say  that  a  function  space,  which  is  defined  by  conditions  on 
the  wavelet  coefficients,  is  robust  if  this  definition  is  independent  of  the  (smooth 
enough)  wavelet  basis  which  is  chosen. 

Let  us  now  motivate  the  formulation  of  the  definition  of  oscillation  spaces  that 
we  will  adopt  for  p  <  0.  Consider  Condition  (3.10)  for  p  >  0.  As  we  already 
remarked,  applying  this  condition  for  j  =  0,  we  obtain  that,  for  any  j'  >  0,  \c\*  |  <  C 
so  that  /  G  5C00 .  Conversely  the  condition  /  G  B]?00  is  necessary  to  make  sure 
that  the  suprema  in  (3.10)  are  finite,  therefore  we  include  this  global  regularity 
condition  in  the  definition  when  p  <  0. 

Since  we  are  interested  in  local  properties  of  functions,  we  will  first  define  these 
spaces  for  periodic  functions  (or  distributions).  We  will  use  periodized  wavelets 
defined  as  follows.  The  functions 

X)24,7¥i)  [2j(x-i)-k],  j>  0,  ke{  o,...,y-i}d, 

iei,d 

together  with  the  function  £2(t)  =  1  form  an  orthonormal  basis  of  the  space  of 
square  integrable  periodic  functions.  We  keep  the  same  notations  as  before  for 
periodic  wavelets  and  the  corresponding  periodic  wavelet  coefficients,  which  will 
bring  no  confusion. 

Definition  22.  Let  p  <  0,  and  s  G  R.  A  periodic  distribution  /  belongs  to 
Op  if  the  two  following  conditions  hold: 

•  /  belongs  to  5^°° , 

•  for  any  e  >  0,  there  exists  exists  C(e)  and  J(e)  such  that 

Vi  >  J(e)  sp-d)jdPx  <  C{c) 2W 

k 


(4.7) 
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Definition  22  does  not  define  a  vector  space  (the  dpx  are  infinite  for  /  =  0).  Its 
remarkable  property  is  that  it  is  independent  of  the  orthonormal  wavelet  basis  (in 
the  Schwartz  class)  which  is  chosen. 

Theorem  5.  For  any  s£l  and  for  any  p  <  0,  the  space  0/  is  robust. 

Let  c\  denote  a  collection  of  wavelet  coefficients  satisfying  |cx  |  <  C.  Let  (ipx) 
be  a  given  orthonormal  wavelet  basis  in  the  Schwartz  class.  We  denote  by  /  the 
distribution  ^2c\ipx.  If  A  =  (yl(A,  A'))  is  an  infinite  matrix  indexed  by  the  dyadic 
cubes,  g  will  denote  the  distribution  whose  wavelet  coefficients  are 

e\  =  ^2  a!(A, 

A' 

(whenever  these  sums  converge),  and  T  will  denote  the  operator  whose  matrix  in 
the  basis  ipx  is  A;  thus  T  maps  /  onto  g. 

In  order  to  prove  Theorem  5,  some  preliminary  remarks  will  be  useful.  We 
start  by  recalling  a  classical  lemma. 


Lemma  3.  Let  7  >  |a|  and  A  £  M.1 .  There  exists  a  constant  C(d)  (which 
depends  only  on  the  dimension  d)  such  that 

Vj,A  \c’j  k  |  <  C2~a^  =>  VA  \ex\<C{d)\\A\\1C2~aB 

This  lemma  expresses  the  fact  that  operators  whose  matrix  in  a  wavelet  basis 
belongs  to  M1  are  continuous  on  Ca  if  |a|  <  7.  It  is  a  straightforward  consequence 
of  Schur’s  lemma  (see  [87]). 


Definition  23.  Let  e  >  0  and  A'  be  a  dyadic  cube.  The  e-neighbourhood  of 
A'  (which  we  denote  by  Ne( A'))  is  the  set  of  dyadic  cubes  A  such  that 


I  j  ~  f  I  <  o' 

<  22ej>2~j' . 


k 

¥ 


k' 

w 


Let  A'  be  a  fixed  dyadic  cube,  and  let  j  be  such  that  |  j  —  j'\  <  ej' .  The  number 
of  dyadic  cubes  A  of  size  2~J  which  belong  to  the  e-neighbourhood  of  A'  is  bounded 

by 

(4.8)  (2j+222ej'2-j')d. 

Note  also  that,  if  A  does  not  belong  to  the  e-neighbourhood  of  A'  and  if  7  >  1/e2, 
then 


(4.9) 


w27(A,  A') 


2-(27  +  d+l)|j-j'| 

(1  +  2mf(ij')rfist(A,  A'))2A+d+1 


<  w7(A,A')2-2'/£. 


Furthermore,  since  /  £  B q’°°,  |cv|  <  C 1;  thus,  if  7  >  2/e,  it  follows  from  Lemma  3 
that 


£ 

A'£Af«(A) 


(4.10) 


A(A,A')ca, 


<  C(d)  ||  A  ||7  C\2~^f- . 
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Now,  we  can  prove  Theorem  5.  Suppose  that  the  wavelet  coefficients  c\  of  / 
satisfy  (4.7)  and  that  A  is  quasidiagonal.  First,  we  notice  that,  since 

J2^(sp~d)jdpx  <  C{c) ¥j, 

k 


it  follows  that 

so  that,  since  p  <  0, 
(4.11) 


VA,  2 (sP~d)idpx  <  C, 


VA,  dx  >  c(e)1/p2~(—dMj2£j/p. 


Let  A  be  a  given  dyadic  cube.  There  exists  A'  C  A  such  that 

(4-12)  |ca'|  > 

For  each  A  we  pick  one  such  A'  which  we  denote  by  A' (A).  We  will  first  prove  the 
following  lemma. 


Lemma  4.  Suppose  that  /  belongs  to  Osp,  and  let  e  such  that 

1 


(4.13)  0  <  e  <  —  . 

1  '  -  2(|s-d/p|  +  l) 

There  exists  C(e)  >  0  and  J(e,A)  such  that  for  any  j  >  J(e)  and  for  any  A  such 
that  |  A |  =  2-J  , 

d\ 


3A  E  Ne( A'(A))  : 
We  will  denote  this  cube  by  A" (A). 


eA"  > 


4 C{d)  ||  A 


-l 


Proof  of  Lemma  4:  If  Lemma  4  were  wrong,  all  cubes  A"  E  N£(A'(A))  would 
satisfy 

\e\n  |  <  - - - . 

-  AC{d)  ||  T-i  || 

Let 

ei„  =eA»  if  A"  E  tVe(A'(A)) 

=  0  otherwise, 

eXn  =  e\n  —e\,, ,  (c^)  =  A -1  (ej^)  and  (c|)  =  A~1{ex).  Since  A -1  is  almost  diagonal, 
applying  Lemma  3,  we  obtain 

(4.14)  VA"  \cU<^- 

We  have  C'(A)  =  T/  4t_1(A'(A),  Ai)eAl  •  Since  /  E  5o’°°  and  A  is  almost 

Ai£JV«(A'(A)) 

diagonal,  Lemma  3  implies  that 

KI<C(/)CV)  ||  a  || 

(where  2-J1  is  the  width  of  Ai);  (4.10)  implies  that,  if  7  >  2/e, 

\Cl{x)\<C(df  ||  T-1  ||7||  A  ||  C(/)2-J,/e. 

Since  e  is  such  that  l/(2e)  >  |s  —  d/p \  +  1,  using  (4.11)  and  j'  >  j,  we  obtain 

2—  (1/2e)V/  <  2 -(s+1~d/p)j  <  dAC_1/,p; 
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that 


2  .  ,  C(d)*\\A-'UA\\C(f) 


lcA'(A)l  < 


C !/P 


dx  2 


-j'/2e 


Thus,  since  l/(2e)  >  1 
(4-15)  lcA'(A)l  ^ 


.  ,0-f  <?(<02  IU-1  imu  II  c(f)  . 
cvml  <  N  - -  I  d\ . 


Choosing  j  large  enough  makes  the  term  between  parentheses  arbitrarily  small; 
thus,  since  c^'( a)  =  ca'(a)  4"  CA'(A)’  (4-14)  together  with  (4.15)  contradicts  (4.12). 


Let  us  come  back  to  the  proof  of  Theorem  5.  Let  e  satisfying  (4.13)  and 
j  >  J(e,A),  where  J(e,A)  is  as  in  Lemma  4;  finally,  let  l  =  [j(  1  —  3e)].  For  each 
cube  n  of  size  2~l ,  there  exists  A  C  p  such  that  p  also  contains  the  cube  A"  (A) 
suplied  by  Lemma  4.  We  denote  these  two  subcubes  A  =  X(p)  and  A"(A)  =  Ai (p). 
Denote  by  2_jl  the  width  of  A"(A).  Finally,  let 

R{p)  =  sup  \e\n  |  . 

A"CjU 

We  have 

R(a)  ^  |eAi(^)  |  ; 

so  that,  using  Lemma  4, 

YJ^{sp-d)lR(A)P  <  J2^(sp~d)l  |eAlWr 

/i  /X 


<  C{A)Y^^sp~d)l d{ 


X(fi) 


(where  the  sums  are  taken  on  all  dyadic  cubes  p  of  width  2  Jl).  Since 

j  G  +  [6e/]}, 

this  quantity  is  bounded  by 

C(A)  2(sp-d)1  E  E< 


3=1  A 


(where  |A|  =  2 


l-\-[6el] 

=  cia)  E  E  dpx  2  ( sp  ~  2  ( “ sp + ^  d  ~ 1 ) 

j=i  a 

<C{A)26 \~sp+d\d  J22isP~d)jdPx- 

H  j=l  A 

Since  e  can  be  chosen  arbitrarily  small,  Theorem  5  follows. 


Thus 


We  have  proved  that  the  definition  of  the  spaces  0*  for  p  <  0  is  independent 
of  the  wavelet  basis  chosen.  These  spaces  are  defined  for  periodic  functions  (or 
distributions).  If  /  is  defined  on  Md,  we  adopt  the  following  definition. 
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Definition  24.  If  p  <  0,  a  distribution  /  belongs  to  Op  loc  if  /  £  Bjh°°  and  if 
its  wavelet  coefficients  satisfy 

\/K  compact,  Ve  >  0,  2i'sp~d~e^  Ma|p  <  C(K,e). 

\£AjCK 

The  proof  that  this  definition  is  robust  is  the  same  as  in  the  periodic  case;  it 
follows  in  particular  that  the  definition  of  u>f(p)  is  robust  for  p  <  0. 

4.3.  Concavity  of  the  scaling  function.  We  will  now  prove  that  the  func¬ 
tion  u>f,  which  is  defined  on  M*,  can  be  extended  on  M  into  a  concave  function, 
and  we  will  also  establish  some  additional  properties.  These  results  are  expressed 
and  proved  more  easily  using  an  auxilliary  function,  the  oscillation  exponent  Sf(p'), 
defined  as  follows. 


Definition  25.  Let  p'  ^  0.  The  oscillation  exponent  of  /  is  defined  by 

Sf(p')=p'uf  • 

Coming  back  to  the  definition  of  uif(p),  it  is  clear  that 

if  p'  e(0,+oo[  sf(p')  =  sup{s  :  /  E  0\/p,  loc}. 

Definition  26.  The  semi-local  Holder  exponent  of  /  is 

(4.16)  Hf  =  sup  {a  :  /  E  C,“c}. 

Theorem  6.  Let  /  £  L°°;  the  oscillation  exponent  Sf(p')  is  a  concave  function 
on  M+;  therefore,  it  is  right  and  left  differentiable  at  every  point.  Its  right  and  left 
derivatives  belong  to  L°°  and  satisfy 

Vp' E  [0, +oo[,  (s'f)r(p1)  <  d  and  (s'f)i(p')  <  d. 

Furthermore,  Sf(p')  is  increasing  on  M+. 

The  scaling  function  u> j  (p)  can  be  extended  on  M  into  a  concave  and  increasing 
function;  its  derivative  uj'j  is  positive  and  decreasing,  and  satisfies 

(4.17)  lim  ui'f(p)  =  Hf. 

p— >-  +  oo  J 


The  fact  that  Sf(p')  is  increasing  on  1+  is  a  direct  consequence  of  the  third 
statement  of  Proposition  10. 

In  order  to  prove  that  ojj  is  concave  on  M,  we  will  first  prove  that  ojj  is  concave 
on  M+  and  on  M-.  Two  cases  can  then  occur:  Either  u>f(p)  =  —  oo  on  M-  or 
u>f(p)  extends  continuously  at  0,  and  (^/),d(0)  >  (‘^/),(j(0);  in  both  cases,  one  can 
conclude  that  w/  is  concave  on  M.  The  concavity  of  ojj  on  M+  will  follow  from  the 
concavity  of  sj]  in  order  to  prove  it,  we  will  need  the  following  inequality  which 
follows  directly  from  Holder’s  inequality.  If0<p<r<g,  and  if  (c4)  is  a  sequence 
of  real  numbers, 

\  !/r  \  a/p  /_^  \(l-a)/q 

(4.i8)  (£Kr)  <(Ei^ip)  (Ei^i9) 


where  a  is  defined  by 
(4.19) 


1  a  1  —  a 
r  p  q 


Lemma  5.  The  functions  Sf  and  ojj  are  concave  on  M+. 


44 


STEPHANE  JAFFARD 


Proof:  Let  p'  =  1/p ,  q'  —  1/q  and  rf  =  1/r.  By  definition  of  sj  ( p' )  and  sj  ( qf ), 
Ve  >  0,  3C  >  0  such  that,  for  j  large  enough, 

J2dpx<2dj2~ps,ip')j 2ei  and  y  d\  <  2dj2~qSql'q''lj2ej . 


Therefore 


y  ^  <  ^ -Psf(p')j^jar/p  ^2~qSf<'ql'>j/1  a'>r/q  2('d+e'lj . 


By  definition  of  Sf(r'),  there  exists  a  sequence  jn  — >■  oo  such  that 

y  drx>2-rsfir''>jn2(d-£'>jn-, 

|A|  =  2-^ 

it  follows  that 

rsf(r')  +  e  >  arsf(p')  +  (1  —  a)rsf(q')  —  e. 

Since  this  inequality  holds  for  any  e  >  0,  and  since  (4.19)  can  be  rewritten 

r'  =  ap'  +  (1  —  a)q' , 

we  actually  proved  that 

Va  E]0,  1[  if  r'  =  ap'  +  (1  —  a)q' ,  then  Sf  (r1)  >  a.Sf  (p1)  +  (1  —  a)sf  (</); 

which  means  that  Sf(p')  is  concave.  But  a  function  is  concave  if  and  only  if  its 
second  derivative  (in  the  sense  of  distributions)  is  negative.  Since  Sf  (p')  =  p'uif  {p), 
then  Sj(p)  =  therefore  ui'j  is  a  negative  measure,  and  ujj  is  concave  on 

M+. 

Lemma  6.  The  functions  Sf  and  ojf  are  concave  onl". 

Proof:  Let  p,  q,  r  >  0  be  such  that 

1  a  l  —  a 


r  p  q 

applying  (4.18)  to  the  sequences  —  with  exponents  —p,  —q  et  — r,  we  get 

4 

\  ra/p  \r(l-a)/q 

yi4r>(yi4id  (yi4i9' 


and  therefore 


1\  ra  r(l  -  a) 

1  -  J  a  U  T  - '  U)  “ T- 


so  that,  since  —r  >  0, 


4  (  r)  >  «4  Q)  +  (1  -  «)4  Q 


and  Sf  is  concave  on  1  .It  follows  as  before  that  u>f  is  also  concave  on  M  . 


We  will  now  prove  the  concavity  of  ojf  onl. 

First  case:  V A  >  0,  3S1  bounded  and  A n  E  SI  such  that 
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Then,  if  p  <  0, 

\dx\p  >  ‘Z~Apjn 

AeA„nn 

so  that 

— uif(p)  >  —  Ap 

i.e.  cj9(p)  <  Ap.  Since  w#  (p)  =  infcA(p),  it  follows  that 
J  f2  J 

VA  >0,  w/(p)  <  Ap; 

so  that  ojf(p)  =  — oo  Vp  <  0. 

Second  case:  3 A  V A,  dA  >  2_U  We  will  hrst  prove  that  w/  is  continuous 
at  0.  Since  the  function  /  is  uniformly  Holder,  it  follows  that  d\  <  1  for  j  large 
enough.  Assume  that  p  >  0;  for  any  bounded  SI, 

X)  lrfT  <  C 

AeAj-nn 

so  that  ojf(p)  >  0.  On  the  other  hand,  it  follows  from  the  upper  bound  VA,  d\  > 
2_V  that 

2_dj  ^  |dA|p  >  2“Aw'; 

AeAj-nn 

so  that  ojf(p)  <  Ap  and  therefore  lim  u>f(p)  =  0. 

p— >-0  + 

By  the  same  argument,  we  get  that,  if  p  <  0,  ojj  (p)  <  0  and  u>f(p)  >  Ap. 
therefore  w /  can  be  extended  at  0  into  a  continuous  function  satisfying 

“7(0)  =  0. 

P  —  Q 

Let  us  now  prove  that  (wy)^(0)  >  (wp^(0).  Let  p  >  0  and  q  <  0.  Let  r  =  - 

and  r  =  - - r  and  r'  are  positive  and  conjugate.  Let  ui  =  — E3-  }et  S3  be  a 

p  p  —  q 

regular  bounded  domain.  Since  no  d\  vanishes,  applying  Holder’s  inequality,  we 

obtain 

2-d^card(Aj  0  S3)  =2“*'  ^  (dx)w 

AeAj-nn 


<(2-^(Ar)‘"(2-*EO"'')1A 
=  (2-*5>,'),,r  (2-« 

Since  Card(Aj  n  S3)  =  2d;>Vol(Q)  +  o( 2dj),  it  follows  that 

“7fo)  ,  “7(g)  <  Q. 

since  r  =  p/w  and  r'  =  —q/uj,  therefore 

“fjp)  _  “7(g)  <  Q 

p  q  ~  ' 
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If _p  — >■  0  and  q  0,  it  follows  that  0)  —  <  0. 


Let  us  now  prove  the  upper  bound  for  sfj.  By  definition  of  sj(pf ), 

<«/(?')>  feO'1/pl> 

therefore  (using  the  second  embedding  between  oscillation  spaces  stated  in  Propo¬ 
sition  10) 

if  q'<p', 

thus,  by  definition  of  Sf(q'), 

sf  {q')  >  s  —  dp '  +  dq  . 


Therefore 


which  can  be  rewritten 


Sf{q')  >  Sf{p')  -  dp1  +  dq1, 


■sf(p')  -  sf(q') 


<  d; 


p'  —  q' 

hence  the  upper  bound  on  the  right  and  left  derivatives  of  s'j . 
Let  us  now  prove  that  u>f(p)  is  non-decreasing.  Since 


(42o)  "'(p|=s'G)“7'G) 

we  only  have  to  prove  that 

Vp'  >0  Sf  (p1)  —  p' s' f{p')  >  0. 

At  p'  =  0,  the  inequality  holds  because  /  is  uniformly  Holder,  so  that  s/(0)  >  0, 
and,  on  the  other  hand,  s'f  is  bounded  in  a  neighbourhood  of  0;  furthermore,  the 
derivative  of  sj  (p1)  —p' s' )  (p')  is  —p's"{p')  which  is  nonnegative.  Since  w/  is  concave 
on  M,  if  it  is  increasing  on  M+,  it  must  also  be  increasing  on  1".  Therefore,  it  is 
increasing  on  M. 


We  checked  that  uj'j  is  nonnegative.  Since  ui"  is  nonnegative,  it  follows  that  uj'j 
is  nondecreasing;  therefore,  it  has  a  limit  when  p  — >■  +oo.  Since  s' j  is  bounded  in 
a  neighbourhood  of  0,  (4.20)  implies  that 

lim  oj'f(p)  =  lim  s/(p). 

p— )-  +  oo  J  p—y  0 


Let  s/(0)  =  lim  s/(p);  if 

p— >-0 

£2'  is  such  that  £2'  C  £2, 


s  <  Hf,  then,  for  any  bounded  £2,  /  £  Cs(£2)  so  that,  if 

sup  |ca'|  <  C2_SJ  . 

A'eA^nO' 


Therefore 

2(-P-rf)j  dxp  <  oo, 

AeAjnn' 


so  that  Sf(j)  >  s,  and  s/(0)  >  Hf. 

Conversely,  if  s  >  Hf,  3£2  bounded  and  jn  — >■  +oo  such  that 


sup  |ca|  >  Aj n2 

AeA^nn 
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with  Ajn  — >■  oo ;  therefore, 

2  (s'p~d)jn  <  2  (S'p-d~SP)A  j{P 

AeAJrenn' 


which  tends  to  +oo  if  s' p  —  d  —  sp  =  0,  hence  if  s'  =  j  +  s.  Therefore  s/  (^-)  <  s  +  j 
Vs  >  Hf.  Letting  p  — y  +oo,  it  follows  that  s/(0)  <  if/;  hence  (4.20)  holds. 


4.4.  Notes.  The  problem  of  finding  a  correct  definition  for  the  scaling  function 
when  p  is  negative  has  attracted  a  lot  of  attention  among  physicists,  especially  in 
the  context  of  fully  developed  turbulence,  see  for  instance  [7,  25,  90,  97]  where 
possible  extensions  of  the  scaling  function  for  p  <  0  are  proposed.  Up  to  now, 
mathematically,  the  scaling  function  was  not  defined  in  terms  of  oscillation  spaces 
Op,  but  in  terms  of  Besov  spaces  5*’°°.  Recall  that  /  belongs  to  Bp’°°  if  its  wavelet 
coefficients  satisfy 

(4.21)  sup2(sp-dW  £  M’’<c 

J  |  A  |  =  2  —  J 

using  this  definition  for  p  <  0  is  clearly  absurd;  indeed  the  quantity  (4.21)  is 
totally  unstable  for  p  <  0  because  one  wavelet  coefficient  can  take  an  arbitrarily 
small  value  ‘accidentally’.  This  explains  why  (to  our  knowledge)  mathematical 
extensions  of  Besov  spaces  to  negative  values  of  p  have  never  been  proposed.  One 
way  to  eliminate  this  source  of  instability  is  to  replace  in  (4.21)  the  single  value 
|cA|  by  a  supremum  of  the  |cA'|  where  A'  is  close  to  A.  This  is  consistent  with  the 
purpose  of  deriving  spectra  of  singularities;  indeed,  a  very  small  wavelet  coefficient 
is  not  the  signature  of  a  large  Holder  exponent  if  it  has  a  large  coefficient  in  its 
immediate  vicinity.  On  the  contrary,  a  small  value  of  the  supremum  means  that, 
indeed,  a  whole  set  of  wavelet  coefficients  close  to  each  other  takes  small  values, 
which  is  the  signature  of  a  smooth  zone.  Therefore,  the  idea  of  taking  suprema  is 
a  sensible  way  to  ‘renormalize’  the  divergence  of  (4.21)  by  disregarding  the  wavelet 
coefficients  which  are  small  ‘by  accident’  (which  means  here  that  there  is  a  large 
wavelet  coefficient  in  the  neighbourhood).  This  is  the  starting  point  of  the  use 
of  the  Wavelet  Maxima  Method  in  this  context,  see  [7,  80,  97]:  The  pointwise 
value  of  the  continuous  wavelet  transform  is  replaced  by  a  supremum  on  all  lines 
of  maxima  ending  at  the  point  considered.  The  wavelet  maxima  method  may  be 
mathematically  unstable,  see  [50];  however,  we  note  that,  in  the  discrete  setting 
supplied  by  orthonormal  wavelet  bases,  it  would  amount  to  replacing  in  (4.21)  the 
coefficient  c\  by  d\  =  supA/cA  |cA'|,  hence  to  replacing  Besov  spaces  by  oscillation 
spaces  in  the  definition  of  the  scaling  function.  Therefore,  it  is  not  surprising  that, 
though  a  scaling  function  based  on  Besov  spaces  has  no  natural  extension  for  p  <  0, 
a  scaling  function  based  on  oscillation  spaces  has  a  natural  and  robust  extension. 

The  idea  of  proving  that  a  criterion  does  not  depend  on  a  particular  wavelet 
basis  by  proving  the  invariance  of  this  criterion  under  the  application  of  an  element 
of  Mi  was  introduced  by  Y.  Meyer  in  [87].  The  notion  of  robustness  was  introduced 
in  [56,  57]  where  it  is  systematically  used  as  a  tool  to  determine  the  largest  possible 
information  that  can  be  derived  from  wavelet  coefficients  and  is  independent  of  the 
wavelet  basis. 
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5.  Prevalent  results  in  multifractal  analysis 

Many  results  of  multifractality  have  been  obtained  for  specific  functions  or 
random  processes.  In  order  to  obtain  ‘generic’  results  of  multifractality,  we  will 
first  discuss  what  is  meant  by  generic. 

5.1.  Genericity  in  infinite  dimensional  spaces.  In  Md,  two  notions  of 
genericity  are  widely  used:  The  notion  of  Lebesgue  almost  everywhere,  and  Baire’s 
notion  of  quasisure.  These  two  notions  share  some  basic  properties  which  are  a 
natural  requirement  for  any  notion  of  genericity  used  in  infinite  dimensional  spaces: 

(1)  Invariance  with  respect  to  dilations  and  translations, 

(2)  Stability  with  respect  to  countable  intersection, 

(3)  Stability  with  respect  to  inclusion  (if  A  is  a  generic  set,  A  C  B  =y  B  is 
a  generic  set), 

(4)  A  generic  set  is  dense. 

Besides  Lebesgue’s  and  Baire’s  notions  of  genericity,  there  exists  many  other 
examples  in  Md;  let  us  mention  the  one-parameter  family  defined  as  follows: 

Let  8  E  [0,d).  A  set  A  is  Hausdorff-8  generic  if  the  complement  of  A  has 
Hausdorff  dimension  less  than  8. 

Let  u  be  a  notion  of  genericity  defined  on  a  vector  space  E  (which,  to  avoid 
trivialities,  is  supposed  to  be  of  dimension  at  least  1);  we  assume  that  the  four 
requirements  listed  above  are  fulfilled,  and  we  denote  by  Gw  the  collection  of  all 
subsets  which  are  generic  for  u>. 

The  collection  of  genericities  has  the  structure  of  a  net:  uj  is  said  to  be  stronger 
than  ui'  if  Gw  C  Gu>  (this  means  that  a  set  which  is  generic  for  u>  is  also  generic  for 
uj').  Equipped  with  this  partial  ordering,  any  family  of  genericities  {w}uEo  has  a 
supremum:  The  genericity  y  defined  by 

AeGfi  if  A  e  Gw  Vw  e  o. 

The  two  strongest  notions  of  genericity  clearly  are: 

•  The  trivial  genericity,  where  the  only  generic  set  is  E  itself. 

•  The  countable  genericity,  where  the  generic  sets  are  the  complements  of 
the  countable  sets. 

The  Hausdorff  genericities  can  also  be  defined  on  infinite  dimensional  spaces.  In 
infinite  dimension,  another  slightly  weaker  notion  is  supplied  by  the  compact  gener¬ 
icity:  The  generic  sets  are  the  complements  of  the  countable  unions  of  compact  sets. 
This  notion  is  clearly  weaker  than  the  previous  ones,  but  stronger  than  Baire  gener¬ 
icity,  and  it  is  usually  too  strong  to  be  useful  in  practice. 

In  finite  dimension,  Baire  genericity  cannot  be  compared  with  either  Lebesgue 
genericity  or  the  Hausdorff  genericities  for  the  partial  order  defined  above;  indeed, 
there  exist  residual  sets  in  the  sense  of  Baire  which  have  dimension  0;  in  space 
dimension  1,  it  is  the  case  for  instance  of  the  set 

A  =  lim  sup  ( —  —  e~q —  e~q 
p£h,q>  0  \9  9 

The  importance  of  proving  a  generic  result  lies  in  the  fact  that  it  implies  a  kind 
of  uniqueness:  If  a  property  V  is  generic,  another  property  Q  incompatible  with  V 
cannot  be  generic  too. 
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5.2.  Prevalence.  Lebesgue  genericity  plays  a  special  role  in  finite  dimension 
because  the  generic  sets  are  obtained  as  the  complements  of  the  measure-zero  sets, 
for  a  “canonical”  measure  which  is  both  cr-hnite  and  shift-invariant:  The  Lebesgue 
measure.  Therefore,  a  natural  question  is  to  wonder  if  such  a  measure  also  exists  in 
an  infinite  dimensional  Banach  spaces.  Unfortunately,  the  answer  is  negative:  There 
does  not  exist  a  cr-hnite  translation-invariant  measure  in  any  infinite  dimensional 
normed  space.  Indeed  a  ball  B  of  radius  1  contains  an  infinite  number  of  disjoint 
balls  of  radius  1/4;  therefore,  if  the  measure  of  B  is  finite,  the  measures  of  the 
smaller  balls,  which  are  all  the  same,  must  vanish.  (The  property  of  the  unit  ball 
in  infinite  dimensional  that  we  used  is  well  known;  indeed,  it  is  the  first  step  of  the 
standard  proof  that  it  is  not  compact:  One  construct  an  infinite  sequence  of  points 
xn  in  this  ball  such  that  Vn  yl  m,  ||  xn  —  xm  ||>  1/2.) 

However,  this  remark  does  not  kill  any  hope  for  an  infinite-dimensional  exten¬ 
sion  of  the  notion  of  translation-invariant  “Lebesgue  measure  zero”;  indeed,  let  us 
consider  the  following  characterization  of  the  Lebesgue  measure. 

Lemma  7.  In  M.d,  a  Borel  set  S  has  Lebesgue  measure  zero  if  and  only  if  there 
exists  a  compactly  supported  probability  measure  p  such  that 

(5.1)  \fxeMd  p(x  +  S)  =  0. 

The  proof  of  this  lemma  is  straightforward:  First,  if  E  has  Lebesgue  measure 
0,  one  can  use  for  p  the  Lebesgue  measure  on  the  unit  cube,  which  clearly  satisfies 

(5.1) .  Conversely,  suppose  that  S  is  a  Borel  set,  and  that  (5.1)  holds.  Then 

0  =  J  p(S  —  x)dx  =  J  meas(S  —  y)dp(y)  =  p(\SLd)meas(S) 

(recall  that  meas  denotes  the  Lebesgue  measure),  so  that  S  has  Lebesgue  measure  0. 

The  characterization  of  the  sets  of  vanishing  Lebesgue  measure  supplied  by 
Lemma  7  does  not  refer  explicitely  to  the  Lebesgue  measure;  therefore  it  can  be 
turned  into  a  definition  in  infinite  dimension  spaces;  the  sets  thus  defined  are  called 
Haar-null,  and  the  notion  of  genericity  that  it  yields  is  called  prevalence.  From 
now  on,  we  will  only  use  this  notion  of  genericity.  The  following  definition  was 
introduced  by  J.  Christenssen  [26]. 

Definition  27.  Let  E  be  a  complete  metric  space.  A  Borel  set  A  C  E  is 
Haar-null  if  there  exists  a  compactly  supported  probability  measure  p  on  E  such 
that 

(5.2)  Viel'1  p{x  +  A)  =  0. 

A  subset  4  of  F  is  Haar-null  if  it  is  included  in  a  Haar-null  Borel  set.  The 
complement  of  a  Haar-null  set  is  called  a  prevalent  set. 

If  (5.2)  holds,  the  measure  p  is  said  to  be  transverse  to  A. 

Remark:  Recall  that  a  probability  measure  is  called  tight  if  Ve  >  0  there 
exists  a  compact  set  Ke  such  that 

p(I<e)  >  1  -  e. 

If  A  is  a  separable  complete  metric  space,  every  probability  measure  on  E  is  tight, 
see  [93].  It  follows  that,  if  E  is  separable,  we  can  drop  the  assumption  that  p  is 


50 


STEPHANE  JAFFARD 


compactly  supported  in  Definition  27. 

The  basic  properties  of  prevalence  and  several  applications  are  detailed  in 
[26,  46]  and  we  refer  to  them  for  additional  information.  The  following  results  show 
that  prevalence  indeed  supplies  a  natural  generalisation  of  the  notions  of  Lebesgue- 
measure  zero  and  almost  everywhere  in  infinite  dimensional  complete  metric  sep¬ 
arable  spaces.  Furthermore,  this  proposition  implies  that  prevalence  satisfies  the 
assumptions  for  genericity  that  we  listed  above. 

Proposition  12.  (1)  If  S  is  Haar-null,  Mx  6  E,  x  +  S  is  Haar-null. 

(2)  If  A  is  Haar-null,  VA  ^  0,  A  A  is  Haar-null. 

(3)  If  E  is  finite-dimensional,  S  is  Haar-null  4=>-  meas(S)  =  0. 

(4)  A  is  prevalent  =y  A  is  dense. 

(5)  A  countable  intersection  of  prevalent  sets  is  a  prevalent  set. 

The  first  and  second  points  assert  that  the  notion  of  Haar-null  (and  therefore 
of  prevalence)  is  translation  and  dilation  invariant,  which  immediately  follows  from 
Definition  27.  The  third  point  is  just  a  restatement  of  Lemma  7.  The  fourth  point 
is  also  straightforward;  indeed,  we  can  clearly  assume  that  p  has  a  support  included 
in  a  ball  of  arbitrarily  small  radius  (by  considering  a  finite  covering  of  the  support 
of  fi  by  balls  Bi  of  radius  e  and  then  correctly  renormalizing  one  of  the  /iIb,  which 
has  a  non-zero  mass);  therefore,  a  Haar-null  set  cannot  contain  a  ball.  The  stability 
by  countable  intersection  is  proved  in  [26,  46].  However,  we  will  give  a  simple  proof 
this  result  in  the  case  where  E  is  separable. 

We  will  say  that  almost  every  element  of  E  satisfies  a  property  V  if  the  set 
of  elements  satisfying  V  is  a  prevalent  set;  equivalently,  we  will  say  that  V  holds 
almost  everywhere  in  E. 

Two  techniques  for  proving  that  a  set  is  Haar-null  are  used:  the  probe  tech¬ 
nique,  and  the  stochastic  process  technique.  Let  us  describe  them. 

If  one  uses  for  transverse  measure  the  Lebesgue  measure  on  the  unit  ball  of  a 
finite  dimensional  subset  V,  Condition  (5.2)  becomes 

Vr  6  F,  (x  +  V)  fl  A  is  of  Lebesgue  measure  zero. 

In  this  case  V  is  called  a  probe  for  the  complement  of  A. 

As  a  simple  illustration,  let  us  prove  the  following  result  which  shows  that 
prevalence  is  a  weaker  notion  than  compact  genericity. 

Proposition  13.  If  E  is  an  infinite  dimensional  complete  metric  space,  then 
any  compact  set  K  is  Haar-null. 

Proof  of  Proposition  13:  Let  /  :lxAx  K  — >  E  be  defined  by 

f(a,  x,  y)  =  a(x  -  y). 

Let  us  first  check  that  Im(f)  E.  Indeed, 

Im(f)  =  {Jf([-N,N]x  I<  x  I<); 

thus  Im(f)  is  a  countable  union  of  compact  sets,  which  are  thus  closed  and  of 
empty  interior;  therefore  Im(f)  is  a  set  of  first  category  which,  by  Baire’s  theorem, 
differs  from  E.  Thus,  let  v  ^  and  let  V  be  the  one-dimensional  vector  space 
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generated  by  v.  We  will  check  that  Mx  £  E,  (x  +  V)  fl  K  has  at  most  one  point. 
Indeed,  if  this  intersection  contained  two  distinct  points  y\  and  1/2,  then 


so  that  v  = 


V 1  ~  2/2 
Ai  —  A2 


3Ai,  A2  :  2/1  =  Xiv  +  x  and  y2  =  \2v  +  x; 

£  Imf;  hence  a  contradiction  follows.  Therefore 
meas((x  +  V)  fl  K)  =  0. 


Recall  that  a  random  variable  X  taking  values  in  E  is,  by  definition,  a  measur¬ 
able  function  X  defined  on  a  probablility  space  IP)  and  taking  values  in  E. 

Such  a  random  variable  defines  a  probability  on  E  by  the  formula  (if  A  is  a  Borel 
set  of  E : 

P(A)  =  P({w  £  Q  :  X{u)  £  A})  =  P(A  £  A). 

Using  this  probability  measure  in  the  definition  of  a  Haar-null  set,  we  see  that,  in 
order  to  prove  that  a  set  A  £  E  is  Haar-null,  it  is  suffucient  to  check  that 


VfeE,F(A+f)  =  0, 


i.e.  that 

V/£U,P(A£  A  +  /)  =  0. 


Suppose  ow  that  V  is  a  property  satisfied  by  some  elements  of  E,  and  that  E  is 
a  space  of  functions  defined  on  (in  that  case,  X  is  called  a  stochastic  process 
if  d  =  1  and  a  random  held  if  d  >  2).  In  order  to  prove  that  the  set  of  functions 
satisfying  V  is  Haar-null,  it  is  sufficient  to  exhibit  a  stochastic  process  X  such  that 


V/  £  E,  a.s.  f  +  X  does  not  satisfy  V. 


We  will  refer  to  this  way  of  proving  that  a  set  is  Haar-null  as  the  stochastic  process 
technique. 

As  an  example  of  application,  let  us  check  which  regularity  results  can  be  ob¬ 
tained  using  for  transverse  measure  the  Wiener  measure  on  the  space  of  continuous 
functions.  The  stochastic  process  A  associated  with  the  Wiener  measure  by  the 
argument  described  above  is  the  Brownian  motion  B.  It  is  easy  to  check  that,  for 
every  continuous  function  /,  then 

Ve  >0,  a.  s.  /  +  Bx  is  nowhere  C1^2+£ . 

Applying  the  line  of  argument  developed  above,  it  follows  that  Va  >  1/2,  for 
almost  every  /  £  C([0,1]),  /  is  nowhere  C1/2+e.  The  use  of  the  Wiener  measure 
cannot  yield  a  better  result;  however,  if  instead  of  the  Brownian  motion,  one  uses 
a  fractional  Brownian  motion  B13  with  0  <  j3  <  1/2,  one  can  show  that,  for  every 
continuous  /,  then 

Ve  >  0,  a.  s.  f  +  B13  is  nowhere  Cl3+e . 


This  yields  an  alternative  proof  of  the  following  result  of  B.  Hunt  concerning  the 
prevalent  Holder  regularity  of  functions  in  Cs ,  see  [45] : 

Proposition  14.  If  s  >  0,  almost  every  function  in  Cs(x 0)  satisfies 

Vt  £  Md,  hf(x)  =  s. 
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Let  us  now  prove  the  stability  of  prevalent  sets  by  countable  intersection  in  the 
case  where  the  space  E  is  separable.  Of  course,  it  suffices  to  prove  the  stability  of 
Haar-null  sets  by  countable  union.  Let  (Afn)ngK  be  a  collection  of  Haar-null  sets. 
To  each  of  them  is  associated  a  random  variable  Xn  defined  on  a  probability  space 
Qn,  taking  values  in  E,  and  such  that 


(5.3) 


V/  E  E  a.s.  /  +  Xn  (/  An. 


Let  SI  =  n“i  ,  endowed  with  the  product  measure  (which  amounts  to  consider 
independent  copies  of  the  random  variables  Xn).  Since,  almost  surely,  Xn  E  E, 
there  exist  Fn  C  S!n  and  N (n)  such  that 

VOef„,  F{dist{Xn,  0)  >  AT(n))  < 

nz 


Let  Y  =  (n-t-JV(n))3  '  By  the  Borel-Cantelli  lemma,  this  series  converges 

almost  surely  in  E.  Let  now  F  E  E  and  m  E  LJ  be  given.  Then 


/  +  Y  —  gm  + 


N(m)): 


-Xm  where  gm  =  f  + 


y  — 

(n 

n  -£m  v 


■N(n)): 


:Xn 


But  gm  belongs  to  E\  since  (5.3)  holds  for  |-m+jv1|.m^3  Xm  instead  of  Xm,  applying 
(5.3)  to  f  +  Y,  we  obtain  that,  for  almost  every  u>  E  Yin^m  ^n,  f  +  Y  Am.  Thus, 
on  the  product  space  SI  =  S 1m  x  n 

V/  E  E,  Vra  ^  N,  a.s.  f  +  Y^Am, 


which  means  that  the  union  of  the  Am  is  a  Haar-null  set. 


We  will  prove  the  following  result,  which  describes  the  prevalent  Holder  regu¬ 
larity  of  functions  of  Op  when  s  >  d/p  (which  coincides  which  the  space  B/’°°  in 
this  case). 

Theorem  7.  Let  s  >  d/p\  then 

•  the  Holder  exponent  of  almost  every  function  /  of  the  space  Osp  takes 
values  in  [s  —  d/p,  s]  and 

(5.4)  VLf  E  [s  —  d/p,  s] ,  df(H)  =  Hp  —  sp  +  d ; 

furthermore,  for  almost  every  x,  hj(x)  =  s; 

•  let  xo  be  an  arbitrary  given  point  in  Md;  then,  for  almost  every  function 
in  Osp ,  hj(x 0)  =  s  —  d/p. 

The  idea  of  the  proof  of  Theorem  7  is  to  find  appropriate  probes  for  Oscillation 
spaces.  We  will  explicitly  construct  bases  of  these  probes  by  defining  their  wavelet 
coefficients. 

5.3.  The  prevalent  spectrum.  In  this  section,  we  prove  the  first  point  of 
Theorem  7.  Thanks  to  the  Sobolev  embeddings,  if  s  —  d/p  >  0,  functions  in  Osp 
cannot  have  Holder  exponents  less  than  s  —  d/p. 

Let  l  E  N  and  M  =  2dl .  We  will  construct  an  M-dimensional  probe  in  Osp.  The 
M  generators  of  this  space  are  defined  by  their  wavelet  coefficients.  Let  j  >  1  and 
k  E  {0,  ...,  2dj  —  1};  K  and  J  <  j  are  defined  by 

k  K 


where  K  E  Ed  —  27Ld. 
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Each  dyadic  cube  A  is  now  split  into  M  subcubes  with  side  2  J  1 .  For  each  index 
i  E  {1,  M}  we  choose  a  different  subcube  i(A).  Let 

(5.5)  a  =  —. 

P 

The  probe  is  spanned  by  M  functions  gi  with  the  following  wavelet  coefficients 


(5.6) 


VA,  d) 


(A) 

dl 

ax, 


=  0  if  A'  is  not  of  the  form  i( A). 


Since  Osp  coincides  with  the  Besov  spaces  Bsp’°°  (because  s  —  j  >  0),  it  follows  that 
the  functions  {<7s}s={i,...,M}  belong  to  0*. 

Let  a  >  1;  a  point  xo  E  is  said  to  be  a-approximable  (by  dyadics)  if  there 

exists  a  sequence  (  Jn,  Kn)  E  x  Ed  such  that 


(5.7) 


x0 


E  n 
2 


< 


1 

2“J- 


(clearly,  we  can  assume  that  Kn  E  ^  —  2Zd). 


Lemma  8.  If  xo  is  a-approximable,  then  there  exists  a  sequence  (rf®.  k)  in  the 
cone  of  influence  of  width  2l  above  xo  such  that 

(5.8)  Vi  \d)tk\>c{M)ra2~H\ 

with  H{=  H (a))  =  s-  |  +  A. 

Proof:  Suppose  that  xo  is  a-approximable  and  let  A  be  the  dyadic  cube  such 
that  and  j  =  [aJ„]  (where  Jn  and  Kn  are  given  by  (5.7)).  The  wavelet 

coefficient  indexed  by  i(A)  has  size 

d\{x)  =  j~a 2(7-*)J'2-pj“  >  c{M)j~a 2(?-*)J'2-£j'. 

Let  /  be  an  arbitrary  function  in  Op  with  wavelet  coefficients  c\.  Consider  the 
affine  subspace  of  dimension  M  composed  of  the  functions 


M 


fp  =  f+'52?9i 


i  =  1 


where  /3  =  (/31 ,  ...,/3M).  Let  xq  E  and  7  >  0.  If  fp  is  (C,  7)  smooth  at  xq  then, 
inside  the  cone  of  width  2l  above  xq, 


(5.9) 


M 

c\  +  'y  ( /3l  d 
8  =  1 


<  C  c(M) 2 


-13 


Denote  by  E“  the  set  of  points  xq  such  that 


3k  : 


x0  - 


2i 


< 


2«j ' 


(Note  that  xq  is  a-approximable  if  xq  E  Ea  =  lim  sup  Ej  ) . 
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The  set  EJ  is  the  union  of  2dj  cubes  of  width  2-2  "U  Suppose  that  x  and  y 
are  two  points  in  the  same  cube  and  suppose  furthermore  that  fp  is  (C,  7)  smooth 
at  x  and  fp  is  (C,  7)  smooth  at  y.  Then  \x  —  y\  <  2  •  2_"J  and  Vi  =  1,  M 

M 

C*(A')  +  I>T(ao  <Cc(M) 2-™' 

(5.10)  <; 

M 

c*-(a')+E^’4a')  <Cc(M) 2-™' 

8  =  1 

for  any  dyadic  cubes  A'  at  scale  j'  inside  the  cone  of  width  2  above  x  and  y.  But, 
since  \x  —  y\  <  2  •  2~aj ,  we  can  find  such  a  A'  satisfying  j'  =  [aj].  Using  Lemma  8 
and  (5.8)  it  follows  from  (5.10)  that 

\\(]  -  p\\  <2Cc(M)2-^~h^  (j')a(-=  A(j')) 

(where  ||/3||  =  sup  \/3i\).  Therefore  the  set  of  [3  satisfying 

i  =  l,...,M 

3x  G  Ej  such  that  fp  is  (C,  a)  smooth  at  x 

is  included  in  the  union  of  2dj  balls  with  radii  A(j').  It  follows  that  the  Lebesgue 
measure  of  the  M-uples  [3  satisfying 

3x  G  Ea  such  that  fp  is  (C,  a)  smooth  at  x 

is  bounded  by 

OO 

(5.11)  ^(Uc(M)[aj]a)M  2-^~H)Mlaj]2dj 

j=J 

(where  J  can  be  chosen  arbitrary  large);  we  can  choose  M  large  enough  so  that 
d  —  (7  —  H)Ma  <  0;  thus,  when  J  tends  to  00,  (5.11)  goes  to  0,  so  that  it  vanishes. 
Therefore,  the  set  of  M-uples  (3  =  (j3i ,  ...,  (3m)  such  that  f  +  ^2,l3%g%  is  (U,  7)  smooth 
at  a  point  in  Ea  has  measure  zero.  Since  it  is  true  for  all  C  >  0,  the  set  of  [3  such 
that 

3x  G  Ea  :  f  +  Y l  /?¥  is  C1 
also  has  measure  zero.  Therefore 

Va  >1,  V7  >  H(a),  a.s.  in  Osp,  V*  G  Ea  hj{x)  <  7. 

Taking  7„  — >  H (a)  (with  7„  >  H (a))  it  follows  by  countable  intersection,  that 

(5.12)  Ma  >  0,  a.s.  Mx  G  Ea  hj(x)  <  H(a). 

Therefore,  if  a„  is  a  dense  sequence  in  (1, 00),  using  the  same  argument,  one  obtains 
that 

a.s.  in  Op,  Vn,  Vx  G  EUn ,  hf(x)  <  H(an)  (V). 

Let  /  be  a  function  such  that  (V)  holds.  Let  a  be  fixed  and  a  subsequence 

of  a„  such  that  is  non  decreasing  and  tends  to  a.  The  subsets  are 

decreasing  and  their  intersection  (:=  Ea)  contains  Ea.  Therefore  any  x  G  Ea 
satisfies  hf(x)  <  H(a)  and  thus 

(5.13)  MxeEa,  hj  (x)  <  H (a). 


WAVELET  TECHNIQUES  IN  MULTIFRACTAL  ANALYSIS 


55 


But  (  see  [52])  there  exists  a  measure  ma  supported  on  Ea  such  that  any  set  E 
of  dimension  less  than  d/a  satisfies  ma(E)  =  0,  and  ma(Ea)  >  0.  Moreover,  if 
Gh  =  :  hf(x)  <  H},  then 


(5.14)  V/  G  0*  dim# (Gh)  <  Hp  —  sp  +  d. 

In  particular,  if  Fa  denotes  the  set  of  points  where  hj  (x)  <  H  (a),  Fa  can  be  written 
as  a  countable  union  of  sets  with  dimension  less  than  d/a  (it  is  the  case  because 
H(a)  =  s  —  |  so  that  ^  =  pH  (a)  —  sp  +  d).  It  follows  that  ma(Fa)  =  0  and 

ma(Ea  —  Fa)  =  0;  but  Ea  —  Fa  is  a  set  of  point  where  the  Holder  exponent  is 
exactly  H(a).  Thus 


mh  g 


df(H)  =  Hp  —  sp  +  d, 


and  (5.4)  holds  on  a  a  prevalent  set. 

Moreover  E 1  =  [0,  l]d,  so  that  we  can  take  mi  equal  to  the  Lebesgue  measure, 
and  (5.13)  yields,  if  a  =  1, 

(5.15)  a.s.  Mx  G  [0,  l]d,  hj(x)  <  s. 


Furthermore,  as  before,  almost  every  function  /  of  O *  satisfies 

meas({x  :  hj(x)  <  s})  =  0; 


so  that 


(5.16)  a.s.  for  almost  every  x  in  [0,  V\d  hj(x)  =  s. 

Results  (5.15)  and  (5.16)  are  not  specific  to  the  unit  cube,  but  they  also  hold  for  any 
cube.  By  countable  intersection,  it  follows  that,  almost  surely,  Mx  G  Md,  hf(x)  <  s 
and,  almost  surely,  a.e.,  hf(x)  =  s.  Therefore  the  first  point  of  Theorem  7  holds. 


5.4.  Regularity  at  a  fixed  point.  In  order  to  prove  the  second  point  of 
Theorem  7,  we  will  use  the  following  regularity  condition,  which  is  a  slight  variant 
of  the  usual  pointwise  Holder  condition: 

/  is  Q7og)- i(*o)  if  there  exists  C  >  0,  d  >  0  and  a  polynomial  P  of  degree  at 
most  [a]  such  that 

if  \x-x0\ <  6,  \f(x)  -  P(x  -  x0)\ <  C- — ^  - p-. 

log  (l/|a?  —  a?0 1) 

A  straightforward  adaptation  of  the  proof  of  the  the  direct  part  of  Theorem  1  shows 
that,  if  /  G  QL,)-i  (*o),  then  there  exists  C  >  0  such  that 

(5.17)  Mj  >  2,  dj(x 0)  <  C—~. 

Let  [3  =  s  —  d/p ;  we  define  a  function  g  by  its  wavelet  coefficients  dj tk  as  follows: 
For  each  j,  only  one  dj  &  does  not  vanish,  and  the  corresponding  hj  is  such  that 
|2Ja?o  —  |  <  2,  in  which  case 

djtkj  =  2_U 

Clearly,  g  G  O/.  The  probe  used  is  the  one-dimensional  subspace  spanned  by 
g.  Let  /  be  an  arbitrary  function  in  Osp.  Let  us  assume  that  there  exist  Ai  and 
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A2  such  that  /  +  Ai g  and  /  +  A2 g  both  belong  to  i(*o)-  Using  (5.17),  the 

wavelet  coefficients  of  /  +  Ai g  and  of  /  +  A2 g  inside  the  cone  \2J xo  —  k\  <  2  satisfy 

1  \  i  I  ,  2-/3j 

\Cj  tk  H-  ^  1  dj  ( k  |  ^  C  -j 

\Cj  ,k  +  2  dj  (  &  |  ^  C  "y  . 

In  particular 

2 -ft/ 

I  Ai  —  A2 1 1  j  I  <  (c  +  c'  ) — — . 

It  follows  that  Vn  >  0,  |Ai  —  A2 1  <  1/ j  so  that  Ai  =  A2.  Thus,  each  line  /  +  Xg  has 
at  most  one  A  such  that  /  +  Xg  belongs  to  C^log^_1  (xo).  Hence,  the  second  point  of 
Theorem  7  holds. 

5.5.  Notes.  The  notion  of  Haar-null  sets  was  introduced  as  early  as  1972  by 
J.  Christenssen  in  [26];  however  this  very  promising  notion  did  not  receive  the  full 
attention  it  deserved  untill  it  was  reintroduced  in  1992  by  B.  Hunt,  T.  Sauer  and  J. 
Yorke,  see  [46].  Several  alternative  notions  of  genericity  have  been  introduced  since 
then,  see  for  instance  [68]  and  the  references  mentioned  there.  The  results  exposed 
in  this  section  are  part  of  a  joint  work  with  Aurelia  Fraysse,  see  [36].  Up  to  now  it 
was  commonly  believed  among  mathematicians  and  physicists  that  multifractality 
was  the  signature  of  very  peculiar  properties  of  the  function  considered  (such  as 
self-similarity  for  instance) ;  therefore  Theorem  7  reverses  the  common  point  of  view 
in  this  held.  The  reader  should  pay  special  attention  to  the  position  of  “almost 
every”  in  the  statements  of  Theorem  7;  indeed  Fubini’s  theorem  does  not  apply 
in  prevalence:  If  one  considers  a  “generic”  function  in  O/,  its  Holder  exponent  is 
almost  everywhere  s,  but  when  a  point  xo  is  fixed,  the  regularity  at  xo  of  almost 
every  function  /  will  be  as  bad  as  possible,  i.e.  s  —  d/p.  Note  that,  in  the  previous 
case,  this  exponent  was  the  one  taken  the  most  exceptionally  (on  a  set  of  dimension 
zero).  The  first  and  second  points  of  Theorem  7  coincide  with  the  Baire-type  results 
of  [55].  However,  here  is  an  example  where  prevalent  and  quasi-sure  results  differ: 
Let  p  <  00;  in  O/, 

•  almost  every  function  satisfies  almost  everywhere  hf(x)  =  s, 

•  quasi  every  function  satisfies  quasi  everywhere  hf(x)  =  s  —  d/p, 

which  follows  immediately  from  juxtaposing  the  results  of  Theorem  7  with  the 
corresponding  results  for  Baire  genericity  in  [55].  We  saw  in  Section  3.2  that,  for 
every  function  of  O * , 

df(H)  <  Hp  —  sp  +  d. 

Therefore,  Theorem  7  shows  that  a  ‘generic’  function  in  O *  is  as  irregular  as  pos¬ 
sible. 

We  obtained  a  generic  result  of  multifractality  in  a  fixed  oscillation  space,  and 
when  s  >  d/p.  It  is  interesting  to  notice  that  the  information  supplied  by  the 
scaling  function  for  p  >  0  can  be  rewritten  as  stating  that  /  does  indeed  belong  to 
an  intersection  of  oscillation  spaces.  Indeed,  it  follows  from  Definition  15  that,  if 
p  >  0,  the  scaling  function  can  be  given  the  following  interpretation:  for  a  given  p, 
if  u>f  (p)  =  u> ,  then 


Ve>0,  /eoj- 


and  / 1 
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Conversely,  if  io(p)  is  a  scaling  function,  let  us  define  the  function  space  Ow  by 

(5.18)  =  5U°  n  (  n  o^pJJpA  . 

\p,  e>0  J 

The  space  Op  is  clearly  a  Banach  space  if  p  >  1;  it  is  a  complete  metric  space  if 
0  <  p  <  1  for  the  metric  defined  by  (3.11).  Because  of  the  concavity  of  u>,  the 
intersection  in  (5.18)  can  be  written  as  a  countable  intersection;  therefore  Ou  is  a 
complete  metric  vector  space.  It  would  be  interesting  to  determine  if  there  is  also 
a  generic  result  of  multifractality  in  Ou ,  and  if  the  multifractal  formalism  holds  in 
this  space. 


The  reader  should  note  that  prevalent  results  cannot  be  obtained  in  the  spaces 
Op  when  p  is  negative  because,  in  this  case,  it  is  no  more  a  vector  space;  this  is 
clear  if  one  adopts  for  Op  the  definition  supplied  by  Definition  22,  since  we  noticed 
that  the  function  0  does  not  satisfy  this  definition.  One  could  however  pick  the 
following  alternative  definition  when  p  is  negative:  /  £  Osp  if  /  £  B^°°  and  if 


(5.19) 


Ve  >  0  3C(e),  J{c) 


Vj>  J(e), 


sp-d+e)j 


1/p 

<C(e). 


This  definition  is  pertinent  only  if  s  >  0;  indeed,  if  /  £  5^°° ,  d\  <  C,  so  that 


i  Ip 

<  C' 2(-s+e/p^ 


which  is  therefore  bounded  if  s  <  0.  This  condition  coincides  with  Definition  22 
when  p  is  positive,  and  therefore  also  supplies  a  natural  extension  of  oscillation 
spaces  when  p  is  negative.  This  time,  the  function  0  belongs  to  the  space,  and  this 
condition  clearly  defines  a  cone;  however,  it  is  not  a  convex  cone  (and,  a  fortiori,  not 
a  vector  space).  Indeed,  let  us  assume  that  s  >  0  and  let  a  be  such  that  0  <  a  <  s. 
In  the  periodic  one-dimensional  setting,  we  define  / 1  by  its  wavelet  coefficients  as 

cjtk  =0  ifO<V<i 


and  f2  by 


_  2-«j 


if  L  <  JL  <  1- 
11  2  ^  ^  ’ 


cj,k 


2~ai  if  0  <  £  <  \ 


=  0  if  \  <  V  <  1. 

The  functions  / 1  and  f2  have  vanishing  d\  as  soon  as  j  >  2,  and  therefore  belong 
to  all  spaces  Osp  for  p  <  0,  but  /i  +  f2  does  not  belong  to  Op  since  a  <  s. 


Even  if  no  prevalent  result  could  make  sense  in  this  setting,  one  could  wonder 
if  there  can  be  a  generic  result  in  some  other  sense;  we  will  check  that  it  is  not 
the  case  for  any  of  the  notions  we  considered;  indeed  all  these  notions  are  based, 
at  least,  on  the  use  of  a  topology  on  the  space  considered;  and  we  will  show  that 
(5.19)  (put  together  with  the  condition  /  £  does  not  allow  to  define  a  finer 

topology  than  the  topology  of  5^°° .  Indeed,  any  given  neighbourhood  of  0  for  this 
finer  topology  would  contain  a  set  of  the  form  Be  P\  Ev,  where  Be  denotes  the  ball 
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centered  at  0  in  B^p° ,  i.e.  is  the  set  of  functions  satisfying  VA ,d\  <  e;  and  Ev 
denotes  the  set  of  functions  satisfying 


Vj  >2 


<  r] 


(here,  we  assume  again  that  we  are  in  the  periodic,  one-dimensional  setting).  Let 
V  be  an  open  neighbourhood  of  0  for  this  finer  topology.  Let 

r  fell 

Br  =  <  /  G  Br  such  that  cj  k  =  0  if  >  —  >  . 


Clearly,  B\/n  C  B\/n  l~l  E o  so  that  it  is  included  in  V  for  n  large  enough,  and,  since 
V  is  open,  there  exist  e'  and  r)'  such  that,  for  any  n  large  enough,  B\/n  +  Bei  fl  Evi 
will  be  included  in  V  and  will  contain  0.  Let 

-  f  k  1 

Br  =  <  /  G  Br  such  that  cj  k  =  0  if  — r-  <  — 

I  L 

For  the  same  reasons,  for  m  large  enough,  B\/m  is  included  in  Bfj  (~\Evi]  therefore, 
for  n  and  m  large  enough,  B\/n  +  B i/m  C  V.  But,  clearly,  B\/n  =  B\/n  +  B i/n  (by 
splitting  any  function  /  as  the  sum  of  a  function  whose  wavelet  coefficients  vanish 
for  >  \  and  a  function  whose  wavelet  coefficients  vanish  for  ^  <  A).  Thus,  for 
n  large  enough  B\/n  is  included  in  V.  Therefore,  the  finer  topology  actually  is  not 
stronger  than  the  B C°°  topology. 
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